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Abstract. We derive a class of energy preserving boundary conditions for 
incompressible Newtonian flows and prove local-in-time well-posedness of the 
resulting initial boundary value problems, i.e. the Navier-Stokes equations 
complemented by one of the derived boundary conditions, in an L p -setting 
in domains f2 C R™, which are either bounded or unbounded with almost flat 
boundary of class C 3— . The results are based on maximal regularity properties 
of the underlying linearisations, which are also established in the above setting. 



Introduction 

We consider the flow of an incompressible Newtonian fluid with constant density 
p > and constant viscosity p > in a domain Q C R™, which is governed by the 
incompressible Navier-Stokes equations 



( N )So B 



where u denotes the velocity field, S — 2pD — pi is the (total) stress tensor with 
pressure p and D — i(Vu+VM T ) the rate of deformation tensor, / denotes the mass 
specific density of the driving force, e. g. gravity, and uo is the initial velocity. The 
domain is assumed to be a half space, a bent half space, or a bounded domain with 
boundary T := d£l of class C 3_ . The flow is subject to a local boundary condition 
prescribed by the linear operator £>, which may depend on the velocity field u and 
the pressure p. If, for a specific flow, the ratio between the fluid viscosity p and its 
density p is sufficiently large, then the inertia term div(pw ® u) may be neglected 
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and the flow is essentially governed by the incompressible Stokes equations 

pdt.u — divS* = pf in (0, a) x SI, 

divu = in (0, a) x O, 

1 >f ' u ° B{u, p) = on (0, a) x flfi, 

u(0) = uq in O, 

which form a linear system of equations. 



We investigate the solvability of (N)^'^ B (locally in time) and of (S)y ^ 8 



under 

diverse boundary conditions and in the L p -setting with p > 1 sufficiently large in the 
non-linear Navier-Stokes case. For standard boundary conditions at impermeable 
walls like the no-slip, the free-slip or the more general Navier boundary condition, 
but also for the Neumann boundary condition, the (local-in-time) well-posedness of 
these systems is well-known. Let us note that even for generalized Newtonian flows, 
the corresponding Stokes system has the property of L p -maximal regularity as was 
recently shown in [3J. In particular, the Navier-Stokes system for a generalized 
Newtonian fluid with one of the standard boundary conditions above is locally-in- 
time well-posed. 

On the other hand, in computational fluid dynamics the problem of formulation 
of appropriate boundary conditions especially at outflow boundaries appears. For 
this purpose, there is a multitude of so-called artificial boundary conditions (ABCs) 
in use, cf. |19j . but a rational derivation of those is often missing. Even more impor- 
tantly, not much is rigorously known about strong solvability of the corresponding 
Stokes of Navier-Stokes system. In fact, some of these ABCs are actually known to 
lead to ill-posed problems. 

The first aim of the present paper therefore is to give a rational derivation of a 
large class of non-standard boundary conditions, containing several different ABCs, 
by introducing the concept of energy preserving boundary conditions. Secondly, we 
will establish an L p -theory of the Stokes system under (standard and) non-standard 
boundary conditions and apply this to the corresponding Navier-Stokes systems, 
obtaining the local-in-time well-posedness of the underlying incompressible New- 
tonian flows. Let us note that L p -maximal regularity of the Stokes system under 
the diverse boundary conditions does not immediately follow from known results 
as, e.g., provided in [9l [10]. The reasons for this are that (i) the known relevant 
maximal regularity theory only applies to parabolic systems, but not to the Stokes 
system and, (ii) the known version of the Lopatinskii-Shapiro condition requires all 
components of the system of boundary conditions to be of the same order. Note 
that, while mixed order parabolic systems are a field of active research, cf. [TT], no 
complete generalization of the Lopatinskii-Shapiro condition for mixed order initial 
boundary value problems seems yet to be available. Therefore, a main part of the 
present paper is to develop the L p -maximal regularity theory of the Stokes prob- 
lem with various boundary conditions. This will be done in such generality that 
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the above mentioned standard boundary conditions are also covered, so that the 
present paper provides a rather complete picture. 

To establish L p -maximal regularity we employ a localization procedure. As a 
starting point, we consider the prototype geometry for initial boundary value prob- 
lems, a halfspace. Using L p -maximal regularity for parabolic systems with diver- 
gence type boundary conditions, which are rigorously derived in appendices, we are 
able to derive the boundary symbols of the halfspace problems and employ state 
of the art methods such as the "H 00 -calculus to obtain L p -maximal regularity. The 
case of a bounded domain with sufficiently smooth boundary is then reduced to 
finitely many (bent) halfspace problems, where the overall approach generalizes the 
well-known techniques for parabolic problems to the Stokes equations. This way, 
we obtain a new generic localization procedure, which is applicable for all boundary 
conditions under consideration. However, the above procedure relies on the ability 
to reduce a fully inhomogeneous Stokes system to the case of a Stokes flow driven 
only by the boundary conditions. Therefore, we also introduce a new splitting 
scheme, which is applicable for all boundary conditions and all geometries under 
consideration. The overall approach may also be applied to study incompressible 
Newtonian flows under dynamic energy preserving boundary conditions of relax- 
ation type or under further non energy preserving boundary conditions; both cases 
are not contained in the present paper. 

The paper is organized as follows. We start with a rational derivation of a large 
class of non-standard boundary conditions in ScctionQ] The main results are stated 
in Section f2J which also includes an overview on known results from the literature. 
The above announced splitting scheme is presented in Section[3J which also contains 
a brief overview of the (weak) Dirichlet and Neumann problems for the Laplacian in 
several geometries, since these are the main tools to achieve the desired reduction 
of the Stokes system. The proof of our main theorem concerning the L p -maximal 
regularity of the Stokes equations is carried out in Sections 0J [5] and [B] for the case 
of a halfspace, a bent halfspace and a bounded smooth domain, respectively. The 
non-linear problem is treated in Section [7] The paper closes with two appendices 
dealing with parabolic systems with divergence type boundary conditions, which 
play a key role for the splitting scheme and the treatment of the halfspace problems. 

1. A Class of Energy Preserving Boundary Conditions 

The derivation of suitable boundary conditions for flow problems is not obvious 
and depends on the physics which is to be modeled there. However, one set of 
equations strictly applies at the boundary T = dfl. These are the balance equations 
at r, also termed transmission conditions. For the relevant quantities mass and 
momentum, these transmission conditions read 



(1) 



\p(u — Ur)] ■ v — 7T on (0, a) x T, 



(2) 



Ipu ® (u — u r ) — Sj ■ v = fr on (0, a) x T, 
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where [</>] denotes the jump of a quantity (j) across the interface if Y is passed in 
the direction opposite to the interface normal v. In these equations u denotes the 
velocity and uy is the interface velocity, which is zero here, since we consider fixed 
domains for which Y is independent of time. As introduced above, S denotes the 
total stress tensor and the right-hand sides in (JlJ and ([2]) model sources and sinks 
due to processes on the boundary, e.g. due to transfer of mass from the bulk to the 
boundary (so-called adsorption). In the situations we have in mind, ?t = but /r 
may be non-zero, motivated by free liquid surfaces with surface tension. 

Now, if r is a physical boundary, by which we mean that it separates the do- 
main SI from a different external bulk phase f2 cx t, then physically sound boundary 
conditions are sometimes available from knowledge about the velocity field in J7 e xt- 
Usually, this will provide one condition and the remaining ones need to be modeled 
based on constitutive assumptions. For instance, if the exterior phase is a rigid 
solid, then M ext = and ([l} implies 

(3) u ■ v = on (0, a) x Y, 

expressing the fact that the "wall" Y = dfl is impermeable. 

To obtain the required additional boundary conditions, consideration of the ki- 
netic energy is helpful. The incompressible flows under consideration are governed 

by 

(4) pd t u + div(pu <g) u — S) — pf, divw = in (0, a) x O. 
Hence the total kinetic energy satisfies 

( 5 ) 5PM 2 dz + J \p\u\ 2 {u-v)da + 2p J \D\ 2 dx - J u-Svda = 0, 
a r n r 

where v : Y — > R™ denotes the outer unit normal field of SI. According to this 
energy balance, the rate of change of total kinetic energy plus the loss of kinetic 
energy due to internal friction equals the power 

(6) 7r NS = j (u- Sv - \p\u\ 2 (u ■ ufj da, 

r 

which changes the total amount of kinetic energy of the system via the bound- 
ary. A boundary condition which causes this contribution via the boundary to 
vanish may therefore be considered as an energy preserving boundary condition for 
incompressible Newtonian flows. 

Consider first the case when Y is impermeable, i.e. <j3j> holds. Then 



ttns = J u ■ Svdtj 
r 

and ((HJ implies 

u ■ Sv — Pyu ■ PySv. 
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where Pr = I — v <g) v. Note that at fixed x 6 T, the map Pr(x) is the orthogonal 
projection onto the plane tangent to Y at x. Evidently, we obtain 7Tns = if the 
no-slip boundary condition 

(7) P T u = on (0, a) x Y 

holds. Together with ([3]), this yields the homogeneous Dirichlet boundary condition 

(8) u = on (0, a) x T 

which has already been introduced by Stokes; see [58] . Let us note in passing that 
the no-slip condition at fixed walls is much under debate recently, especially in the 
context of flow through micro-channels; cf. [5J. Its inhomogeneous version, i.e. 

(9) u = tiD on (0, a) x Y, 

is used in numerical simulations to model inflow boundaries, where the velocity 
profile needs to be known within sufficient accuracy. 

Another obvious way to have 7Tns = 0, assuming ([3]) to hold, is the free-slip ( or, 
perfect-slip) boundary condition 

(10) PvSv = on (0, a) x Y. 
The complete boundary condition then reads 

(11) u ■ v = and 2^iP T Dv = on (0, a) x Y. 

Note that at points where T is locally planar, a simplification is possible. Indeed, 
in this case the tangential derivatives of the velocity have no normal component 
due to ([3]), hence 2/j.PrDi' reduces to /iPrd^u and, hence, the complete boundary 
condition becomes 

(12) u ■ v = and [iPyd v u = 0. 
In the applied literature this is often written as 

u- v = and d v U\\ = 0, 

which this has to be understood in the right sense: at a fixed point x G Y it means 
d y {u-T) = at x for every fixed tangential vector r JL v{x). Since d v {u-r) — d v u-T 
for fixed r, this is indeed equivalent. Note that for non-planar Y, the condition (fT2")) 
differs from (jlll) by additional curvature-related terms. 

The free-slip and the no-slip conditions are the two extreme cases a = and 
a — > oo, respectively, of the more general Navier boundary condition 

(13) P T u + a P r Sv = with a > on (0, a) x Y. 

This condition, which is due to Navier [3S], implies 7Tns < 0, but 7Tns does not 
need to vanish. Note that, mathematically, the Navier condition is a lower order 
perturbation of the free-slip condition. Therefore, it will not play a role later on. 
There is a more general Navier-type partial-slip condition, which has recently been 
shown in [S] to be in some sense the most general boundary condition which is 
possible for flows of incompressible Newtonian fluids at impermeable walls. 
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Another standard boundary condition is motivated by ([2]) as follows. If f2 C xt 
is a gas phase with negligible gas viscosity, then a reasonable simplification of the 
momentum transmission condition is 

-Sv = p cxt v + f r on (0, a) x T 

with p cx t the external pressure. The associated homogeneous boundary condition 
is the homogeneous Neumann boundary condition 

(14) Sv = on (0, a) x T, 

which models a free surface with external pressure set to zero and vanishing surface 
tension (/r = 0). Under this condition, energy need not be preserved for the 
Navier-Stokes system, but for the Stokes equations. Indeed, we have ng = 0, where 



(15) 7rs = / u-Svda. 

r 

To obtain further boundary conditions, which are energy preserving for the 
Stokes system, we split both u and Sv into normal and tangential parts and obtain 

(16) 7T S = / {{u ■ v) (Sv ■v)+P v u- PrSvJ da. 

r 

Hence another admissible combination is ([7]) together with Sv • v = 0, i.e. 

(17) P r u = and 2(id„u ■ u - p = on (0, a) x T. 
In the applied literature it is usually written in the form 

(18) P T u = and 2fid u (u ■ v) - p = on (0, a) x T 

and is employed as an outflow boundary condition; cf. the remark behind (fT2j). 
Often, the factor 2 in front of the viscous term is omitted; cf. [TH] and note that the 
factor 2 does not appear if one derives a kinetic energy balance backwards, starting 
with the Navier-Stokes system in which div S has been replaced by [iAu — Vp, 
employing already divu = 0. 

For incompressible flow and planar outflow boundary, the condition Pru = 
together with divu = implies d v (u ■ v) = for sufficiently regular solutions. This 
way one obtains another outflow boundary condition, namely the pressure condition 

(19) Pru = and p = Po on (0, a) x T. 

The latter two boundary conditions are examples of so-called artificial boundary 
conditions (ABCs), which are imposed at artificial domain boundaries, being inside 
the flow domain. 

To motivate further boundary conditions, which are employed as ABCs in the 
numerical literature, we first need the following simple observation. Due to the 
incompressibility condition, it follows that 

div D = div R = ^ Au, div S — div T = fiAu — Vp, 
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where R = |(Vu — Vit T ) denotes the rate of rotation tensor (or, spin tensor) and 
T = 2/iR — pi is the antisymmetric counterpart of the stress tensor S. Hence, the 
alternative form 



— J ^p\u\ 2 dx + J ^p\u\ 2 (u ■ v) da + 2p J \R\ 2 dx — J u-Tvda 



= 0, 



of the kinetic energy balance is also available. Consequently, the power 7Tns can also 
be expressed as 

(20) 7Tns = / (u ■ Tv — 5/c|it| 2 (u ■ v)j da. 



If this contribution vanishes, the total kinetic energy will also be monotonically 
decreasing and we hence also consider boundary conditions which imply zero power 
due to expression (I2D1) . 

Starting again with the case of an impermeable wall T, i.e. assuming ([3]) to hold, 
a further energy preserving boundary condition evidently is 

(21) u-v = and 2pRv = on (0, a) x T. 

In K 3 , a direct computation shows that this is in turn equivalent to rot u ■ r = for 
any r _L v. This leads to the vorticity boundary condition 

u ■ v = and Prrc-t u = on (0, a) x T. 

For the special case of the Stokes system, additional energy preserving boundary 
conditions can be read off. The power input via the boundary then is 

(22) 7T S = J u ■ Tv der, 

r 

and u ■ Tv can be decomposed according to 

(23) vr s = / ( P r u ■ P r Tv - (u ■ v) p) da; 



r 

note that Tv -v = —p. Hence, replacing ([3]) by the complementary equation Tv-v — 
leads to the Neumann- type boundary condition Tv = 0. In the split form it reads 
as 

(24) -p = and 2pRv = on (0, a) x T, 

and it may be considered as the homogeneous version of a vorticity-pressure bound- 
ary condition, related to a free boundary. 

There is one more admissible combination, namely 

P r u = and Tv ■ v = on (0, a) x T. 

This leads to the homogeneous version of the outflow boundary condition (fl9)) . i.e. 

(25) P r u = and - p = on (0, a) x T; 

but note that this time it also appears in this form for general curved boundaries. 
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Finally, we will consider the well-posedness under two further boundary condi- 
tions which are not energy preserving, even for Stokes flow, but appear naturally 
as combinations of partial boundary conditions given above. These conditions arc 
also employed as ABCs in the numerical literature and they read as 

(26) 2[iP T Dv = and - p = on (0, a) x T, 

which corresponds to the prescription of the tangential part of the normal defor- 
mation rate and the external pressure, and 

(27) 2(iRv = and 2p, d v u ■ v - p = on (0, a) x T, 

which corresponds to the prescription of the tangential part of the vorticity and 
the normal part of the normal stress. 

For a better overview of the boundary conditions which are rigorously analyzed 
in this paper, we summarize them below. As energy preserving boundary conditions 
for the Navier-Stokes system we have the conditions (|8|), (fTTT) and (|2~Tj) . which read 

(Bla) Pru = and u ■ v = on (0, a) x T, 

which equals the no-slip condition [u] = 0, 

(Bib) 2 t iP T Dv = and u ■ v = on (0, a) x T, 

which is known as the perfect slip condition, and, finally, 

(Blc) I^lRv = and u ■ v = on (0, a) x T, 

Additional energy preserving boundary conditions for the Stokes system are the 
conditions |14]), (JT7J) , ((24]) and ([25]), which read 

(B2a) Ptu = and 2[i d v u ■ v — p = on (0, a) x T, 

which corresponds to the prescription of tangential velocities and the normal com- 
ponent of normal stress, 

(B2b) 2 f iP r Dis = and 2fi d v u ■ u - p = on (0, a) x T, 

which equals the Neumann condition Sv = 0, 

(B2c) P r u = and -p = on (0, a) x T, 

which corresponds to the prescription of tangential velocities and the external pres- 
sure, and, 

(B2d) 2iiRv = and - p = on (0, a) x T. 

Finally, we also consider the two boundary conditions (|26[) and (f27|). which read 

(B3a) 2[iP T Dv = and - p = on (0, a) x T, 

which corresponds to the prescription of the tangential part of the normal defor- 
mation rate and the external pressure, and 

(B3b) 2fiRv = and 2(1 d v u ■ u - p = on (0, a) x T, 
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In the remainder of this paper we show that all of these conditions lead to locally- 
in-time well-posed Navier-Stokes systems in the appropriate L p -setting. 

2. Main Results 

Our analysis of incompressible Newtonian flows subject to one of the boundary 
conditions (B) is based on i p -maximal regularity of the underlying linear system, 
i. e. the incompressible Stokes equations 

pdtu — /j,Au + Vp = pf, div u — g in (0, a) x fi, 

( S )/f,f « P) = h on (°> °) x dfl > 

it(0) = Mo in O 

with fully inhomogeneous data /, g, h and uo- Here B denotes the linear operator, 
which realizes one of the discussed boundary conditions (B). 
We will focus on the cases where fi is the half-space 

M™ := { (x, y) G R" -1 x M : y > } , 

a bent half-space 

K ■= { (*> V) e R"" 1 x M : y > w(x) } , 

with a sufficiently smooth and flat function uj : R n_1 — > M, or a bounded domain 
with sufficiently smooth boundary. In all cases we require the boundary T — d£l of 
51 to be of class C 3 ~ and we assume I < p < oo. 

To establish maximal regularity in an L p -setting, we employ the natural solution 
spaces 

u G X» := ^((0, a), L„(n, R n )) n L„((0, a), # p 2 (f>, R")), 
where if* denotes the Bessel potential space of order s > 0, and 

peX p (a) :=L p ((0, a), H^Q)), 

where 

Hp(Q) := {4>&V(Q) : V(f> G L p (Ct, M") } 

denotes the homogeneous Bessel potential space of order one, which becomes a 
semi-normed space with 

Wm p (n) ■■= l|V^|| M n,Kn), 4>eHfa). 

Note, that Hp(Sl) is a dense subspace of Hp(Q) for all domains under consideration. 
As usual, T>'(Q) denotes the space of distributions on £1 Note, that the regularity 
assumptions on the domain f2 imply the embedding 

ffi(Q) ^ { G L Ptloc {VL) : 4> G SJ'CSl open and bounded } 
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to be valid for all 1 < p < oo, cf. [36., Chapitre 2, Theoreme 7.6]. The corresponding 
data spaces are therefore determined as 

/ e Y f (a) := L p ((0, a), L p (Sl, R™)), 

.9 G Y 9 (a) := V^ 1/2 ((0, a), L p (fi))nL p ((0, a), fl^fi)), 

u e Y u := V^ 2/p (r!, R") 

and the regularity class for h depends on the boundary condition. 

Note, that in the L p -setting it is necessary to carefully distinguish between a 
function defined on 17 and its trace on T. To account for this, we denote by [•] the 
trace operator. 

To simplify our notation, we denote the boundary operator as B = B a '^ , where 
the first index a G { — 1, 0, 1 } is used to describe the tangential part of the bound- 
ary condition and the second index j3 G {—1,0, 1} is used to describe the normal 
part of the boundary condition, i. e. a — indicates the prescription of the tan- 
gential velocity Pr [u] , a — 1 indicates the prescription of the tangential part of the 
normal stress 2^,Py[D]v, a — —1 indicates the prescription of the vorticity 2/i[i?]^, 
(0 = indicates the prescription of the normal velocity [u] • v, j3 = 1 indicates the 
prescription of the normal part of the normal stress 2fid u u ■ v — \p], and, finally, 
P = — 1 indicates the prescription of the pressure — [p\. Thus, the parameters a 
and j3 will be used to describe the order of the corresponding part of the boundary 
condition as \a\ resp. \/3\ and to simultaneously fix the particular convex combi- 
nation between Vu and Vu T , which would have to be used to obtain a boundary 
condition of order one based on the kinetic energy balance or its alternative form. 
Hence, we consider the linear operators 

B a ^(u, p) = P r B a 'P(u 7 p) + Q T B a ^{u, p), 

where Qr = I — Pr denotes the projection onto the normal bundle of T, given as 

(Bl) P T B°^(u, p) := P r [u], P T B ±l ' := ^P T [Vu±Vu T ]v 

for & G{-1, 0, +1} and 

Q r B a '°(u, p)-v:= [u] ■ v 

(B2) 

Q r B a < +1 (u, p)-v.= 2fxd v u ■ v - \p], Q v B a ~ x {u, p)-v:= -\p] 

for a G { —1, 0, +1}. Note, that the normal derivative in the L p -setting has to 
be understood as d u = [V • T ]v. Also note, that the boundary operators B a, ° 
with a G {—1, 0, +1} realize the energy preserving boundary conditions (Bl) for 
incompressible Newtonian flows. Moreover, the boundary operators B 01 ^ 1 with 
a G { —1, 0} and B a ' +1 with a G { 0, +1 } realize the additional energy preserv- 
ing boundary conditions (B2) for incompressible Newtonian Stokes flows. Finally, 
the boundary operators B~ 1,+1 and yB +1,_1 realize the non-preserving boundary 
conditions (B3). 
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Now, if u £ X u (a) and p £ X p (ct), we first obtain by trace theory 
P r i3°' /3 (u, p) G T°(a), where 
T°(a) := ^- 1/2p ((0, a), L P (I\ IT)) n L p ((0, a), W^ 1/p (T, IT)) 
for j3 £ { — 1, 0, 1 }, where Tr denotes the tangent bundle of T, 
P r B ±1Jj (u, p) £ T^a), where 
T± x (a) := ^ p 1/2 - 1/2p ((0, a), L p (I\ Tr)) n L p ((0, a), W^ 1/p (T, XT)) 
for p£ {-1,0, 1}, 

Qr# Q '°(u, p) G N£(a), where 

N°(a) := ^- 1/2p ((0, a), L p (I\ AT)) n L p ((0, a), W^ 1/p (I\ AT)) 

for a £ { — 1,0, 1}, where AT denotes the normal bundle of T, and, finally, 

QrB a,±1 (u, p) G Nf(a) := L p ((0, a), I^- 1/p (r, AT)) 

for /3 e { — 1, 1 }. Therefore, the regularity class for the boundary data ft is given 
as 

ft G Y^(a) := „ e L p ((0, a), i p , ioc (r, M™)) : 

i (»7|i/)i/eNg(a) J 

for a, /3 £ { -1, 0, 1}. 

However, some boundary conditions may impose additional regularity properties 
on the pressure trace [p] . Since 

Q r B a ' +1 (u, p) ■ v — 2/i d v u ■ v — [p] 

for a e {—1, 0, 1}, the regularity class of the data may also be chosen according to 
the regularity of [Vw] . Obviously, the regularity of [Vu] is significantly higher than 
the regularity of [p\. Indeed, as will be shown, in case of the boundary conditions 
given by £> Q >+ n for a £ { —1, 0, 1 } every regularity condition in between that two 
may be imposed, i. e. the pressure p belongs to the regularity class 

X P;7 (a) := { Tr £ X p (a) : [w] £ % 7 ((0, a), L p (0)) n L„((0, a), ^"^(r)) } 

with 7 S [0, 1/2 — l/2p], if and only if the boundary data satisfies 

Qrft G N+ 7 (a) := ((0, a), L p T, AT)) n L p ((0, a), % 1 " 1/p (r, AT)). 

In case of the boundary conditions given by B 01 '^ 1 for a £ { —1, 0, 1 }, the normal 
part of the boundary data equals the trace of the pressure. Hence, the pressure p 
belongs to the regularity class X Pi7 (a) with 7 £ [0, 00), if and only if the boundary 
data satisfies Qrh £ 7 (a) := NjJ"*(a). For simplification we will therefore employ 
the notations 

Xp.-oo^Xpfa) and Nf r J ) := Ng(a) 
for /3 e { — 1, 0, 1 } and require 

ft G Yg(a) := ry e L p ((0, a), i p , ioc (r, M™)) : 

[ (#>e< 7 (a) J 
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with 7 = -oo for /3 = 0, with 7 £ {-00} U [0, 1/2 - l/2p] for f3 = 1 or with 
7 G { -00 } U [0, 00) for /3 = -1. 

Another issue concerning the pressure p, which has to be addressed, is the fact, 
that its uniqueness and continuous dependence on the data can only be guaranteed, 
if the pressure trace [p] is uniquely determined via the boundary condition. If 
the pressure does not appear in the boundary condition, uniqueness has to be 
understood as uniqueness up to a constant. To account for this phenomenon, we 
define 

X p ,-oo( a ) : = X p,-oo(a)/R and X^*(a) := X Pj7 (a), 
which ensures a unique pressure p £ X^ (a) , if the boundary condition is given by 
the linear operator B = B a '@ with a, ft £ {—1, 0, 1}. The additional regularity 
parameter 7 has to be chosen according to the constraints discussed above. Note, 
that X° _ oc (a) constitutes a Banach space, whereas X^ ^^(a) is semi-normed via 

However, for 7 > 0, the spaces X^ 7 (a) equipped with their natural norm 

IMIx±i,(a) = max { klx,(a)> llWIU^(a) } ' 77 G X P,7( a ) 

constitute Banach spaces, too. Analogously, the data spaces Yj'^fa) and Y^'^(a) 
with a G { — 1, 0, 1 }, /3 £ { —1, 1 } and 7 > constitute Banach spaces with their 
natural norm 

IMIy";£(o) = max { H P r?7l| TS(a ), WQtvWmI^v) } ' ^ e Kn( a )> 
whereas the spaces ^^.'-00 with a £ { — 1, 0, l},/3e{— 1,1} are semi-normed via 

l^k.'-ooM = maX { l^k.-ooM } ' 77 £ Y V?( a )- 

Hence, continuous dependence of the solution on the data has in some cases to be 
understood w. r. t. semi-norms, regardless of its uniqueness, which will always be 
guaranteed. 

In addition to the above regularity conditions, there are several compatibility 
conditions, which have to be satisfied by the data. First of all, the compatibility 
condition 

(Cl) ff ,„ divu = 3 (0) 

is necessary and a boundary condition given by the linear operator B = B a, @ with 
a, f3 £ { — 1, 0, 1 } implies the compatibility conditions 

Prho] = Prh(0), if a = and p > ~, 
(C2)£ u 2 
' ° /iiMVwo ± Vu~l)v = P r h(0), if a = ±1 and p > 3 



to be necessary for (S)*'^ ' hua to admit a maximal regular solution. 



Last but not least, there is a n additiona l compatibility condition, which stems 

ZaTTB 



from the divergence equation in (S) 



> f,g,h,u 



To reveal it, we set 
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with l/p+ 1/p' = 1 and define the linear functional 

(•,•): Y g (a) x Nl-ocia) — > L p ((0, a), H-\Q)) 
for %j) G Y ff (a) and n G N° h , .^(a) via 



(01 ("0, J7)l 



J Will da- J H>dx, <j>&H l pl (£l). 



An integration by parts yields 

(^|(div«, [«]•!/))= /V^-udx, 4>eH£,(n) 



gJ^(O) 



and we infer 

1(01 (divu, [u] •!/)>!< Ilull^^l^l^i^n), 

as well as 

|(0|<9 t (divu, N -f ))| < l|w|| Xu (a)l0liji,(a): 

Hence, a boundary condition given by the linear operator 23 = i5 Q '' 3 with a, /3 £ 
{ — 1, 0, 1 } implies the compatibility condition 

<9r[«o] = Qr/i(0), if p > |, and 

(g, Qvh) G flj((0, a), | • 

there exists fj€ _oo(a) such that 
Qr[uo]=r){0), if p > § , and 
( ff ,7 7 )e J ff p 1 ((0, a), (#1(0), I • | A i #(n) )'), 



(C3)^, 



if P = 0, 



if/?G{-i, 1} 



to be necessary for (S) j'^ )f UQ to admit a maximal regular solution. 

On the other hand, the above regularity and compatibility conditions are also 
sufficient to construct a unique maximal regular solution to the Stokes equations 
(S)/g 'h,u f° r a ^ 5' ^' u °) G Y" ,,3 (a), where we choose the space Y" ,/9 (a) to 



consist of all 



(/, .9, ft, u ) G Y/(o) x Y fl (o) x Yg(a) x Y„ 



Q 'f . The detailed results will be 

g,h,u 



which satisfy the compatibility conditions (C) 
stated as Theorem [53] below. 

With the above notations at hand, we are able to formulate our main results. 
Concerning the local well-posedness of incompressible Newtonian flows subject to 
one of the boundary conditions (B) we will prove 

Theorem 2.1. Local Well-Posedness, Semi-flow, Energy Inequality. 

Let Q C R" be a half-space, a bent half-space or a bounded domain. Let V = dfl 
be of class C 3 ~ and let n + 2 < p < 00. Let B = B aJi with a, f3 G { -1, 0, 1 } be 
one of the linear operators (B), which realizes one of the boundary conditions (B). 
If /3 = let 7 = -00; if p = 1 let 7 = 1/2 - l/2p; if = -1 Zei 7 > 0. 



14 



DIETER BOTHE, MATTHIAS KOHNE, AND JAN PRUSS 



Then the Navier-Stokes equations (N)°?^' admit a unique local-in-time solution 
(u, p) on a maximal time interval [0, a*(f, Mo)); whenever the data f G Y^(oo) 
and uq G Y u satisfy the compatibility conditions (C)q'q Uq . The solution is in the 
maximal regularity class 



ueX u (a), peXj i7 (a) 



for all < a < a*(/, uq). The maximal existence time is characterized as 
a*(f, uq) < oo =>■ lira u(t) does not exist in W 2 ~ 2 ^ p {Vt). 

t->-a*(/, Ho) 

Moreover, the solution enjoys the following properties: 

(1) If one of the energy preserving boundary conditions for incompressible New- 
tonian flows (Bl) is imposed and f € Y/(oo)nL p ((0, oo), £2^, K ra )), then 
the energy inequality 



— J ip|u| 2 dx < J \pu\f )dx 



n n 
is valid. 

(2) If f = 0, then the solution map uq m- u generates a local semi-flow in 
Z a >? = { v G W*- 2/p (n) : v satisfies (C)^ v } , 
the natural phase space for (N) 1 ^^ ' in the L p -setting. 



Note, that the main statement of Theorem 12.11 is the existence of unique local- 
in-time solutions as can be seen as follows. 

Remark 2.2. The claimed characterization of the maximal existence time is a direct 
consequence of the existence of unique local- in-time solutions. Indeed, ifa*(/, uq) < 
00 and lim t _j. a ,(y UQ ) u(t) =: v G Wp 2 ^ P (Q) would exist, then there would also be 
a unique local-in-time solution starting at t = a*(f, uq) with initial value v, which 
would extend the solution starting at t = with initial value uq beyond its maximal 
existence time. Analogously, the semi-flow property of the solutions to (N)J5°^ Q ' 
is a direct consequence of the existence of unique local-in-time solutions as well. 
Moreover, the claimed energy inequality is a direct consequence of the construction 
of the energy preserving boundary conditions (Bl). Hence, for a complete proof 
of Theorem 12.11 to be established, it is sufficient to prove the existence of unique 
local-in-time solutions with the claimed regularity properties. 

The construction of unique local-in-time solutions will be carried out in Section[7l 



based on the maximal regularity property of (S)™'^ „ in the L p -setting, which is 
provided by 

Theorem 2.3. L p -maximal Regularity, Semi-flow, Energy Inequality. 

Let a > and let Q C R™ be a half-space, a bent half-space or a bounded domain. 
Let r = dtt be of class C 3 ~ and let 1 < p < oo, p ^ §, 3. Let B = B a ^ with 
a, j3 G {—I, 0, 1 } be one of the linear operators (B), which realizes one of the 
boundary conditions (B). 
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If /3 = let 7 = -oo; if f3 = 1 let 7 e { -00 } U [0, 1/2 - l/2p}; if f3 = -1 lei 
7 G { -00 } U [0, 00). 

Then the Stokes equations (S) a f^'^ UQ admit a unique maximal regular solution 

ueX u (a), peX^(a), 

if and only if the data satisfies the regularity conditions 

feY f (a), geY g (a), h € Yg(a), u eY u 

and the compatibility conditions (C)g'^ Ug . The solution map 

Y^(»)^I a (tt)x^(a) 

is continuous and the solution enjoys the following properties: 

(1) If one of the energy preserving boundary conditions for incompressible New- 
tonian Stokes flows (Bl) or (B2) is imposed, and, additionally, f G Yy(a)n 
L p ((0, a), £2(^7 R")); 9 = and h = 0, then the energy inequality 




n n 

is valid. 

(2) If f = 0, g = and /i = 0, i/ien t/ie solution map uq ^ u generates a 
semi-flow in 

1 a > p = { v E W^ 2/P {n) : v satisfies (C)^ v } , 

the natural phase space for (S)q'^' 0uo in the L p -setting. 

Again, the main statement of Theorem 12.31 is the existence of unique maximal 
regular solutions as can be seen as follows. 

Remark 2.4. Once the existence of unique maximal regular solutions is proved, the 
continuity of the solution map is a consequence of the open mapping principle. 
Moreover, the semi-flow property follows from the existence of unique maximal 
regular solutions. Last, but not least, the claimed energy inequality is a direct 
consequence of the construction of the energy preserving boundary conditions (Bl) 
and (B2). Hence, for a complete proof of Theorem 12.31 to be established, it is 
sufficient to prove the existence of unique maximal regular solutions. 

We close this section with some remarks, how the considered boundary condi- 
tions have already been treated in the literature. Concerning strong solutions to the 
Navier-Stokes equations subject to a homogeneous Dirichlet condition, H. Fujita 
and T. Kato as well as P. E. Sobolevskii have established unique local strong 
solutions in an L 2 -setting based on a semigroup approach already in 60's, cf. [T3irHl 
152"] . Later on, Y. GiGA and T. MlYAKAWA as well as F. B. Weissler generalized 
these results to the L p -setting, cf. [T7J [TH1 [HI] ■ The first approaches based on resol- 
vent estimates in an i p -setting are due to V. A. Solonnikov, M. McCracken 
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respectively S. Ukai, cf. [53j [34l [60] . Finally, maximal Lp-regularity was estab- 
lished by V. A. Solonnikov, W. Borchers and T. Miyakawa, respectively 
W. Desch, M. Hieber and J. Pruss, cf. [53J El H2]- A semigroup approach to 
the Stokes and Navier-Stokes equations, which in particular yields maximal L p - 
-Lg-regularity, was developed by T. Kubo and Y. Shibata, cf. [28l [29j [30] . This 
approach was later generalized to Navier respectively Robin boundary conditions 
by Y. Shibata and R. Shimada, cf. [44 . The extremal cases, i.e. the no-slip and 
the perfect slip condition, are also covered and the results apply to a halfspace, 
bent and perturbed halfspaces as well as to bounded and exterior domains. 

The Stokes and Navier-Stokes equations subject to Neumann boundary condi- 
tions have first been considered by V. A. Solonnikov in a series of publications, 
cf. [54[ [55j [56j [57]. Later on, Y. Shibata and S. Shimizu established maximal 
Lp-regularity, cf. [151 HH1 S3 EE1 EH]- In all cases, the motivation are either prob- 
lems with free boundary or two-phase problem with an evolving phase-separating 
interface. As has already been mentioned in the introduction and in Section [I] 
the Stokes and Navier-Stokes equations subject to a Neumann boundary condition 
arise as a model problem in these situations. 

A totally different approach to the Stokes and Navier-Stokes equations in an L p - 
setting was developed in the 90's by G. Grubb and V. A. SOLONNIKOV, cf. [531 

[HdSl [551 [5D1 [HH2. In thiS Seri6S ° f P ublications > tlic Stokes and Navier-Stokes 
equations are transformed into a system of pseudodiffcrcntial evolution equations 

and treated by an abstract pseudodifferential calculus. This way, the authors are 

able to treat the Dirichlet, perfect slip and Neumann boundary conditions and 

establish maximal Lp-regularity in each case. As a remarkable fact, this method 

is applicable for mixed order boundary conditions, which have not been treated in 

the literature before. 

Another popular boundary condition for impermeable walls, which is not men- 
tioned in Section [1] arises from the Navier condition, if we additionally assume the 
boundary F to be perfectly flat, i. e. to coincide with a two-dimensional plane in 
M 3 . In this case, the Navier condition coincides with the Robin Condition. For 
this condition, a complete L p -theory including an %°°-calculus due to J. Saal is 
available by [42] |43] and the monograph [41] . 

A further well-known approach to the Navier-Stokes equations dates back to the 
fundamental works of J. Leray and E. Hopf, cf. [32 [33] [27] , who introduced the 
concept of weak solutions. This way, one may construct global solutions without any 
smallness assumption on the initial datum. However, the question of the uniqueness 
of these solutions still remains open in space dimension n > 3. Therefore, weak 
solutions to the Navier-Stokes equations are still an active field of research in math- 
ematics and the concept is further developed nowadays. Nevertheless, the literature 
focuses on (homogeneous) Dirichlet conditions. Notable exceptions are the articles 
by H. Bellout, J. Neustupa and P. Penel, cf. [H [371 EHl [Ml ED] , who consider 
boundary conditions of Navier type as generalized impermeability boundary condi- 
tions, and the article by F. Boyer and P. Fabrie, cf. [4 , who consider Neumann 
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type boundary conditions. For a more detailed overview of the theory of weak solu- 
tions and its development we refer to the monographs by O. A. Ladyzhenskaya, 
R. Temam and G. P. Galdi, cf. [3T| [59) \15\ [16] . 

Several of the energy preserving boundary conditions (B) are well-known to be 
applicable for artificial outflow boundaries in numerical simulations. The vorticity 
conditions (a = —1) and the pressure conditions (/3 = —1) have first been used in 
[8]. However, a mathematically rigorous analysis of the resulting initial boundary 
value problems is not available up to now. 

3. A Splitting Scheme 

To construct a solution to (S)^'" ' hUa with B = B 01 ' 13 , where a, f3 G { —1, 0, 1 }, it 
will be convenient, to reduce the problem to the spacial case / = 0, g = 0, Prh — 0, 
uq = and 7 = 1/2 — l/2p, if /3 = 1. To achieve this reduction, we decompose 
the desired solution u € X u (a) and p G X^ 7 (a) as u — v + u with v, u G X H (a) 
and p = q+ p with q, p € X£ 7 (a), where v is determined as the strong solution to 
a suitable parabolic problem with data depending on /, g, Prh and uq, and q is 
determined as the weak solution to a suitable elliptic problem with data depending 
on /, g and Qph. If [5 = — 1, we may choose u = and p = 0, i. e. the solution is 
completely determined by the splitting scheme. If (3 G {0, +1 }, the remaining part 
u and p of the solution is determined as the maximal regular solution to (S)q'^'^ , 
where Prh = 0. Of course, h will depend on v and q, but we will ensure 

Qrh G ^((0, a), W^p(T, NT)) n N^-^a), if /? = 0, 

resp. 

Or^N+^^W, if/3 = +l. 

Before we state our main theorem concerning the splitting scheme, we need to 
prepare the suitable framework for the elliptic boundary value problems, which will 
be solved for the pressure. Depending on the boundary condition under considera- 
tion, we will construct solutions q G Hp(Q) to the Dirichlet problem 

(DP) fh -A<7 = div/ mil, [q] = h on d£l 

with data / G L p (Q, R") and h G Wp' 1/p (dfl). Therefore, we want to recall the 
known results concerning this problem for a sufficiently large class of domains and 
allow f2 to be the halfspace R™, a bent halfspace R™, where uj : R n_1 — > R is 
assumed to be of class C 1 , or a bounded domain with boundary of class C 1 . For 
all these domains the trace space Wp 1 ^ p (dfl) is well-defined and there exists a 
bounded linear trace operator 

(28) {^■.H'W^W^idCl). 

Due to the weak regularity assumptions on q and /, problem (DP) has to be 
understood in a weak sense. To achieve the corresponding weak formulation, we 
define 

oflJ(n) := els (C °°(Q), I • |^ {n) ) C ffi(O) 
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as the closure of the space Co°(f2) of compactly supported, smooth functions in 
the semi-normed space Hp(Cl). Then, for all domains under consideration and all 
1 < p < oo the identity 

^(n) = {^j p 1 (fl) : M = o} 

is valid, cf. [51]. With this definition the (weak) Dirichlet problem (DP) j,/, is 
equivalent to its weak formulation 

-(W|V«) = (W|/), <t>e H^(Sl), 

\ )f,h 

[q] = h on dfl. 

Following the argumentation in [51) . this is exactly the right setting to obtain 
maximal regular solutions to (DP) in a weak sense via 

Proposition 3.1. The Weak Dirichlet Problem. 

Let O C M" be the half-space, a bent half-space or a bounded domain with boundary 
r = dfl of class C 1 and let 1 < p < oo. Then the (weak) Dirichlet problem (DP)f t h 
resp. (DP)™ h admits a unique maximal regular solution q G Hp(Q), whenever f 6 

L p (fl, R") and h G Wp~ 1/p (dQ). 

Proof. First note, that the trace operator (|2"8|) is onto. Now, given / G L p (£l, K ra ) 
and h G Wp -1 ^ (dd), we first choose q G Hp(fl) with [q] — h and then solve 

-( | v? - v? ) = ( | / + V5), € (n) 

to obtain q — q E qHHSI). The possibility of obtaining such a solution follows from 
[5"Tl Theorem II. 1.1], which covers the case of a bounded domain, resp. [STj Lemma 
II.2.5], which covers the (bent) halfspace case. 

Finally, if q, q G Hp (SI) are two solutions to (DP)™ then q—q£ a Hp(Sl) solves 

-(V^|Vg-Vg) = 0, 0e 
which implies g = <j again by [51] Theorem II. 1.1] resp. (5TJ Lemma II.2.5]. □ 

A main ingredient of the proof of the above result was the possibility to construct 
extensions in Hp(SY) to traces in Wp 1 (dSl). On the other hand, Proposition 13. II 
also implies this possibility in an even refined sense. 

Proposition 3.2. The Trace Space of Hp (SI). 

Let SI C K™ be the half-space, a bent half-space or a bounded domain with boundary 
r = dSl of class C 1 and let 1 < p < oo. Then the trace operator 

[■}-.H^(n)^w^(dn) 

is onto and there exist a bounded linear extension operator 

e : w}-v*(m) — > Hp-(n), 

which is characterized by 

-AEh = 0inSl, [Eh] = hondSl, h G Wl~ 1,v (dSl). 
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Proof. Given h £ W p 1 ^ p (dCl), define Eh £ H^fl) to be the unique solution to 
the (weak) Dirichlet problem (DP) 0i /,. □ 

Another valuable consequence of Proposition 13.11 is the validity of the Weyl 
decomposition 

(w) Lp(fi, K n ) = s p (n) e v flj(n) 

for all domains under consideration. Here 

S*p(ft) :={/eL P (0, M") : div/ = 0} 

denotes the space of L p -functions, which are divergence-free in the sense of dis- 
tributions. Since both, S p (fl) and Vo-ffp(fi) are closed subspaces of L p (tt 7 R"), 
this direct decomposition is topological and the thereby induced bounded linear 
projection 

W p : L p (n, R") — ► L p (n, M") 

onto S p (Q) along Voijp(^) is called Weyl projection. Given / £ L p (Sl, R n ), we 
have W p f = f— Vq £ S p (fl), where q £ oH p (fl) is obtained as the unique solution 
to(DP)_ /)0 . 

On the other hand, we may construct solutions q £ H p (^l) '■= H p (£l)/M. to the 
Neumann problem 

(NP) M - A 9 = div / in Q, d v q + [f]-u = h on dVL 

with data / £ L p (n, R") and h £ W p 1/p (dSl) := W p 1 ,~ 1/p ' (00)', where 1/p+l/jf = 
1. The formulation of this problem requires a proper definition of the left hand side 
of the boundary condition, since Vq, / £ L p (Ct, R ra ) do not possess traces in the 
sense of Sobolev spaces. However, we may employ a generalized normal trace 

[■l:S p (n)^w p ^(dn), 

which is defined via 

(1>\Mu) :=(v£viH i)£W p l r 1/p '(dn), V £S p (n). 

Note, that for (j) £ H pl (Q) we have E[<j>] — 4> £ oH pl (Q) and, therefore, 

Wp>(VE[<t>] - V<f>) = 0. 

Hence, 

(VE[4>]\v) = {VE[(j)\\W p v) = {W p ^E[4>]\v) 

= (W p ,Vcj>\v) = (Vcj>\W p v) = (Vcj>\v), <P£H£,(n), v£S p (Q) 
and the generalized principle of partial integration 

(MIK) = (V^|«), 06^(0), ue5 p (fi) 

is available. Especially, we have \ [v] v ) = (V(/)|u), for all ip £ W^~ l ^ p (<9f2), 
4> £ Hh(Cl) with [0] = ip and w £ S p (fl), which implies the definition of the 
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generalized normal trace to be independent of the particular choice of the extension 
operator. The Neumann problem (NP) t ^ has to be understood as 

Vq + feS p (n) and [Vq + f]„ = h on dVt 

and is therefore equivalent to its weak formulation 

(NP)/,h -{V<t>\Vq) = {Vct>\f)-{[<t>]\h), 6^,(0). 

The next proposition, which is proved in |50j . shows, that the (weak) Neumann 
problem also admits maximal regular solutions. 

Proposition 3.3. The Weak Neumann Problem. 

Let f2 C R™ be the half-space, a bent half-space or a bounded domain boundary 
r = dfl of class C 1 and let 1 < p < oo. Then the (weak) Neumann problem (NP)f t h 
resp. (NP)J h admits a unique maximal regular solution q e H p (Cl), whenever f e 
L p (Q, W 1 ) and h G Wp 1/p (dQ). □ 

A first interesting consequence of Proposition 13.31 is the pendant of Proposi- 
tion [321 for the generalized normal trace operator. 



Proposition 3.4. The Trace Space of S P (Q,). 

Let ft C R™ be the half-space, a bent half-space or a bounded domain with boundary 
r = dfl of class C l and let 1 < p < oo. Then the trace operator 

[■l:S p (n)^w-^(dn) 

is onto and there exist a bounded linear extension operator 

E u : W- 1/p (dSl) — ► S p (n) 

with 

E u h G VJ3"jJ(n), [E v h} v =h ondfl, he W p ^ p {dn). 

Proof. Given h G W p l ^ p (dCl), define q G H p (Q) to be the unique solution to 
(NP)o,h and set E v h := Vq. □ 

Last, but not least, Proposition 13.31 implies the Helmholtz decomposition 

L p (fl, R") = S p (fl) © VH${Q) 

to be valid for all domains under consideration and all 1 < p < oo. Here 

oS p (n) :=cls R") : div«£ = 0},| • | ip(w ) C L p (0, R") 

denotes the closure of the space of solenoidal, compactly supported, smooth vector 
fields in L p (Q,, R"), which may be characterized as 

s p (n) = {cf>es p (n) : [0]„ = o>. 

Since both, qS p (£1) and V£Tp(fi) are closed subspaces of L p (Q, R"), this direct 
decomposition is topological and the thereby induced bounded linear projection 

Hp : L p (n, K") — -> L p (n, R") 
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onto o$p(0) along ViI p (Q) is called Helmholtz projection. Given / <E L p (ft, R"), 
we have H p f — f — Vg £ oS p (Q), where q £ Hp(Q) is obtained as the unique 
solution to (NP)_/ :0 - 

With the above definitions and results at hand, we may deal with an inhomo- 



geneous divergence condition in (S)i' ff ' ftuo by defining a suitable pressure, whose 



gradient forces the velocity field to develop the requested divergence. 
Proposition 3.5. A Divergence Adjusting Pressure. 

Let a > and let f2 C R™ be the half-space, a bent half-space or a bounded domain 
with boundary T — dfl of class C 3 and let 1 < p < oo . Let 

g £ % 1/2 ((0, a), L p (n)) n Lp((0, a), 
such that there exists 

T] £ W^ 1/2p ((0, a), L p (dn)) n L p ((0, a), w P 2 - 1/p (dr!)) 
with (g, n) £ H^((0, a), (fl"* (fi), I • lijwn))'), 
where l/p+ 1/p' — 1. TTien £/iere exists 

q £ L p ((0, a), offp 1 ^)) 
wit/i — divVq = (pdt — pA)g in X>'(f2). 

Proof. Since H p ,(Q) is a dense subspace of H^,(Q), there exists a unique extension 
F e H*((0, a), oHp 1 ^)) to ( 5 , 77), where # P (^) := H^(Q)'. Now, we define 
q £ L p ((0, a), ffi(O)) be the unique solution to 

(V0|Vg)H0|F), 0efl£,(n), 

which is available thanks to Proposition 13.31 cf. also [SO]. On one hand, we have 
dtVq £ L p ((0, a), L p (Cl, R™)), and, on the other hand, 

(V</>\Vq) = (<t>\F) = ( ( t>\(g,r ] ))=- / 5 dz, ^€C °°(fi) 

Q 

implies Ag = divVg = 5 and hence AWq = VAg = V.g G £ p ((0, a), L p (f2, R™)). 
Finally, we define 9 e L p ((0, a), H p (fl)) via Vg = — (X — W p )(p<9 f - ^A)Vg and 
infer 

-div V<? = div (1 - W p )(pd t - pA)\7q = div (pd t - /iA)Vg 
= (pd t - /iA) div Vg = (pd t - /xA)g in 2?'(fi). 

□ 



Using Proposition 13.51 we may formulate the promised splitting scheme, which 
may serve to eliminate almost all inhomogeneous data for the Stokes equations 



{ a )f,g,h,v 



Theorem 3.6. A Splitting Scheme. 

Let a > and let f2 C R™ 6e £/ie half-space, a bent half-space or a bounded domain. 
Let r = dQ. be of class C 3 ~ and let 1 < p < 00, p ^ §, 3. Let B = B a '^ with 
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a, f3 G {—1, 0, 1 } be one of the linear operators (B), which realizes one of the 
boundary conditions (B). 

If p = let ~f= -co; if P = 1 let 7 e { -00 } U [0, 1/2 - if P = -1 lei 

7 G {-oo}U [0, 00). 

T/ien ifte unique maximal regular solution 

u G X„(a), p e X^ j7 (o) 
to </ie Stokes equations (Sff g 'huo w tth data 

feY f (a), geY g (a), h G Y^(a), u G Y u 

satisfying the compatibility conditions (C) g f tUo may be obtained as follows: 

1. If P — 0, choose rj = h ■ v, define q G Lp((0, a), oHp(Cl)) via Proposition 
based on g and rj and set q = q + K; if P € { — 1, 1 }, choose r\ G N° _ oc (a) 
according to the compatibility condition ( C3) g h Uq , define q £ L p ((0, a), oHp(Cl)) 
via Provosition \3.5\ based on g and r] and set q = q — E(h ■ v). 

2. Choose v G X u (a) £0 6e the unique maximal regular solution to the parabolic 
problem 

pd t v — fxAv — — Vg + pWpf, in (0, a) x Cl, 
P T B a {v) = P T h, [divv] = [g] on (0, a) x dCl, 
v(0) = uq in Cl. 

and define q G X^ (a) by means of Vq = Vq + p(I — W p )f . Then, if P = — 1, 
u = v and p = q constitute the unique maximal regular solution to the Stokes 
equations (SJ/ g ; h uo . 

3. If P G {0, 1 }, £/ie solution to the Stokes equations (S) a f' g ' huo splits as u — v + u 
and p = q + p, where 

ueX u (a), peX^fo) 

denotes the unique maximal regular solution to the Stokes equations (S)q^j*q 
with h — h — B a '@(v, q). In particular Prh — and 

Q r h g ^((o, a), wv 1/p (r, at)) n K,-M, ifP = 0, 

resp. 

Qrh€ti+\ /2 _ l/2p (a), ifP = +l. 

Before we prove this theorem, some remarks seem to be in order. In step 2 we 
used the abbreviation 

P T B a (v) = P r B a ^(v, q), 

since the tangential part of the boundary condition neither depends on the pa- 
rameter P nor on the pressure q. Analogously, we will make frequently use of the 
abbreviation 

QtB^(v, p) = QrB a 'P{v, q), 

since the normal part of the boundary condition does not depend on the parameter 
a. The parabolic problem employed in step 2 will be treated in the appendix. 
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Indeed, there exists a unique maximal regular solution v £ X„(a), thanks to the 
regularity assumptions on the data and the compatibility conditions (Cl) SjU0 and 
(C2)^ Ua . The overall strategy for the proof of Theorem 13.61 in connection with the 
proof of Theorem 12.31 will be as follows. 

Remark 3.7. In the following proof of Theorem 13.61 we will only show, that the 
solution constructed in steps 1 and 2 has the desired regularity properties and solves 
the Stokes equations, where the normal part of the desired boundary condition may 
be violated if (3 £ {0, 1}. In these cases, the first two steps reduce the task to 
prove the existence of unique maximal regular solutions to the Stokes equations as 
claimed in Theorem 12.31 to the special case / = 0, g = 0, Pph = and uq = with 
Qrh £ ffi((0, a), W- 1/P (T, AT)) n K,-M, * = 0, resp. 7=1/2- l/2p, if 
(3 = 1. Once this additional task is accomplished, Theorem 12.31 and Theorem 13.61 
will be completely proved in these cases, cf. also Remark 12.41 On the other hand, if 
P = — 1, the first two steps establish the existence of a maximal regular solution to 
the Stokes equations and it remains to prove its uniqueness to complete the proof 
of Theorem [531 and Theorem 13.61 However, if [u, p) £ X u (a) x X~ ^(a) constitutes 
a maximal regular solution to the Stokes equations (S)q q q q ano - P ~ ~ 1j t nen 
p £ L p ((0, a), oHp(n)) and the momentum balance delivers 

-(V0|v P ) = o, 

which implies p = by Proposition ^. 31 and, hence, u = by uniqueness of the solu- 
tion to the corresponding parabolic problem. Therefore, the proofs of Theorem l2.3l 
and Theorem 13.61 for the case j3 = — 1 will be complete. 

Proof of Theorem Iff. 61 First note, that by construction we always have 
q 6 L p ((0, a), o-Hp 1 ^)): -div Vg = (pd t - ^A)g in P'(fi) 
and therefore 

q £ X^(a), -div Vg = (p9 t - ^A) 5 in D'(fi) 

and g satisfies the boundary condition 

-[g] = /»•!/, if/3e{-l, 1}. 

Hence, q £ X^ (a) enjoys the same property on the boundary. Moreover, v £ X„(a) 
by the maximal regularity property of the parabolic system and 

pd t v - pAv + Vq = pd t v - pAv + Vq + p(l - W p )f = pf in (0, a) x fl, 

which shows, that the momentum balance is valid for v and q. Applying the diver- 
gence in the sense of distributions to the partial differential equation used to obtain 
v and using the properties of q as well as the compatibility condition (Cl) Sl „ , we 
derive 

pdt(divv — g) — /iA(divu — g) = in (0, a) x O, 
[divw — g] = on (0, a) x <9f2, 
(divu-s)(0) = inQ. 
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Hence, div v — g by uniqueness of weak solutions to the diffusion equation with 
Dirichlet boundary conditions. Last, but not least, v satisfies the tangential part 
of the desired boundary condition and the initial condition v(0) = uq. 

Now, if = —1, v and q constitute the unique maximal regular solution to the 
Stokes equations (S)*'- On the other hand, if € {0, 1}, u = u — v and 
p — p — q solve the Stokes equations 

pdtu — pAu + Vp = 0, div u — in (0, a) x ft, 
P r B a (u) = 0, QtB 13 ^, p) = h on (0, a) x dCl, 
u(0) = in fi 

with h = Qr(h- [«]), if = 0, resp. h = -2pd u v -v, if = 1. □ 

4. The Halfspace Case 

This section is devoted to the proof of Theorem 12.31 in the halfspace case, i. e. 
we assume a > 0, ft = K™ and 1 < p < oo with p ^ |, 3. Following Remark 12.41 
it is sufficient to prove existence of a unique maximal regular solution for all data 
satisfying the stated regularity and compatibility conditions. However, following 
Remark |3"771 we may restrict the proof to the case a e { —1, 0, +1 },/?£{ 0, +1 }, 
/ = 0, .g = 0, Pph = and u$ = 0. Moreover, it is sufficient to consider the shifted 
equations with a — oo, since these also deliver maximal regular solutions to the 
Stokes equations on finite time intervals. 

Due to the simple geometry of the domain we may split the spatial variable into 
a tangential part x £ and a normal part y > 0. Moreover, we may split 

the velocity field as u — (v, w) into a tangential part v : [0, oo) x R™ — > M™ -1 
and a normal part w : [0, oo) x R™ — > R to obtain the system of interior partial 
differential equations 

pev + pd t v — pA x v — pdyV + \7 x p = in (0, oo) x R" , 

pew + pd t w — pA x w — pdyW + d y p = in (0, oo) x R™ , 

\/ x -v + d y w = in (0, oo) x R™, 

which have to be complemented by the given initial and boundary conditions. The 
arbitrary parameter e > denotes the shift. Due to v(0) = and w(0) — we 
may employ a Laplace transformation in time and a Fourier transformation in the 
tangential variable to obtain the transformed system 

oj 2 v - nd%v + i£p = Re A > 0, £ e R™ _1 , y > 0, 

oj 2 w - pid 2 y w + d y p = Re A > 0, £ € R"" 1 , y > 0, 

i£ T u + = Re A > 0, £ G R n_1 , y > 0, 

where u, it) and p denote the transformed components of the solution, A G C 
denotes the Laplace co- variable and £ € R n_1 denotes the Fourier co- variable of x. 
Moreover, we used the abbreviations 

A £ :=e + A, u := ^/ p\ e + [i\£,\ 2 . 
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Now, an exponential ansatz leads to the representation 



v(X, f, y) 
w(X, £, y) 
p(x, & y) 



iC ICI 

kA £ 



"T 



0e" l£l!/ 



Re A > 0, 

£ G M™" 1 , 2/ > 0, 



where we have set £ := ^/JIS, and k := p^/Ji and where z u and z w denote the trans- 
formed components of a function z = (z v , z w ) : [0, oo) X R" _1 — > R™, which has 
to be determined via the boundary conditions. We will treat the various boundary 
conditions under consideration separately in the next subsections. Note, that the 
normal component of the boundary datum h is the only non-zero part of the data. 
In the sequel we will denote it by h w = —h- v = h n . Also note, that every boundary 
condition will lead to a linear system 










_z w (X, £)_ 




_h w (X, 0. 



Re A > 0, £ e 



with a family of linear operators B a '*(X, £), which will be shown to uniquely de- 
termine z(X, for all Re A > and £ E R™" 1 . 

The Case a = and /3 = 0. Here we have to treat the plain Dirichlct conditions 



[v] = 


and 


M = 




on 


(0, oo) 


x 9M™ 


and we may assume h w € 


o^((0, 


oo), Wp 1/P (R 




The boundary conditions lead to 














B°'°(A, C) = 


UJ 


-K ' 
ICI . 




a; 



" 

UJ 




1 

UJ 



Now, we have 
1 

1C_ jci 

UJ 

and, hence, 



UJ 

ICI 



UJ j UJ 



*C T 



&u = S 1 - 



iCJ\ JCI 

UJ J UJ 



- 1 h w = 1- 



ICI 



icrx 



This implies 



9„p 



[^P] = ^A £ |C|^=pAjl 



ICI 



= pA £ 



1+ 1 



^ = w(w + |C|)^ 
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and, hence, the desired solution may be obtained by solving the splitting scheme 
pen + pd t u — pAu = —Vp, —Ap = in (0, oo) x K" 
P r B°(u)^0 1 [divu]=0, d v p = T°h w on (0, oo) x dW% 
u(0) = in M" , 
where the bounded linear operator 

T° : ^((0, a), W-VPQRi-l)) R N o^ M _^ L ^ a)> ^-i/p (R »-1)) 

is defined via its Laplace-Fourier symbol 

m (A, f) =w(w+|C|), 

i. e. T° = (pg + p<9 t - pA r ) 1/2 (pe + p<9 t - 2pA r ) 1/2 , e. g. by the "H°°-calculi of the 
involved operators. Indeed, the calculations in Appendix B imply the solution of the 
above splitting scheme to solve the shifted Stokes equations, cf. relation (To o)- On 
the other hand, any solution to the shifted Stokes equations with h w = satisfies 
the equations of the above splitting scheme with a boundary condition 



d v p - 
-i/p 



on (0, oo) x 



for some h p e L p ((Q, a), W p ' (R n )). Employing relation fTo. of once again, we 



infer h p = T [w] = T°h w = 0. Therefore, p = const, and u = 0. 

The case a = ±1 and f3 — 0. In this case the boundary conditions read 

Tp[d y v] — pV x [w] = and [w] = h w on (0, oo) x cM" 

and we may assume h w € H^((0, oo), MV 1/p (R"- 1 ))nN^ .^(oo) by Theorcm[3J> 
The boundary conditions lead to 

±v^ 2 



B ±1 '°(A, 



V£(<C®iO -VmKICI±<ICI) 



Kl 



±^/JiUJ 2 



o u 



it 



& Ki- 
Ki 



Now, we have 



it 



/ ic Ki ± liiii 1 ) 

V CJ CJ UJ u ' 



KKl- 
Ki 

Co' 

f - 



icn ki 



<®»C ICI 

2 



i£Ki ± KKi 

UJ LU LU UJ 



1 ± 



II- 



and, hence, 



{(i-s)^r(i±s)-^ 



= i- 
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This implies 



3r h w = (u 2 ± \C\ 2 )K 



and, hence, the desired solution may be obtained by solving the splitting scheme 
peu + pd t u — fiAu = — Vp, -Ap = in (0, oo) x K™ 
P T B ±l (u) = 0, [divw]=0, d v p = T ±l h w on (0, oo) x dR" 

u(Q) = inM-l, 
where the bounded linear operators 

T* 1 : 0^(0, a), W~ 1 / p (R n ~ 1 )) n N°(oo) — ► L p ((0, a), ^"^(Rn-i)) 
are defined via their Laplace-Fourier symbols 

m±i(A, 0=^ 2 ±IC| 2 , 

i. e. T+ 1 = (pe + - 2^A r ) and T _1 = {pe + pd t ), e. g. by the "H°°-calculi of the 
involved operators. Indeed, the calculations in Appendix B imply the solution of the 



above splitting scheme to solve the shifted Stokes equations, cf. relation (T^o ). On 
the other hand, any solution to the shifted Stokes equations with h w — satisfies 
the equations of the above splitting scheme with a boundary condition 

fin 



d v p = h p on (0, oo) x 
-Vpi 



for some h p e L p ((0, a), W p 1,1 '(R n ~ 1 )). Employing relation (fT^ 



once again, we 



infer h p = T* 1 ^] = T ±1 h w = 0. Therefore, p = const, and u = 0. 



The Case a = and /3 

\v] = and 



1. In this case the boundary conditions read 

- 2p{d v w} + \p] = h w on (0, oo) x dR^, 

and we may assume 7 — 1/2 — l/2p by Theorem 13.61 The boundary conditions 
read 

-K 



to 



kX £ + 2 V ^|C| 2 



OJ 



Now, we have 
1 

UJ 

and, hence, 



1 -i— — I CI 

2 T 2 TP 



2^/JIuj 2 



ifl I , I Kf 

ui 2 "r" 2 to 2 



2 \/ 



1 
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This implies 



[p] = k\ £ z w = p\ e (l 
= pX 



and, hence, the desired solution may be obtained by solving the splitting scheme 
pen + pdtu — pAu = —Vp, —Ap = in (0, oo) x K" 
P r B°(u) = 0, [div«] = 0, \p]=S°h w on (0, oo) x dW^ 
u(0) = in R™, 
where the bounded linear operator 

5° = Nfc,i/ a -i/2p(«>) — > N^ 1/2 _ 1/2p (oo) 

is simply the identity. Obviously, a solution of the above splitting scheme solves 
the shifted Stokes equations. On the other hand, any solution to the shifted Stokes 
equations with h w = satisfies the equations of the above splitting scheme with a 
boundary condition 

[p] = on (0, oo) x dRl 
and we infer p = as well as u = 0. 

The case a = ±1 and f3 — +1. In this case the boundary conditions read 

=F/i[<9yv] — /LtV K [«;] = and — 2p[d y w] + [p] = h w on (0, oo) x <9R" 

and we may assume 7 = 1/2 — l/2p by Theorem 13.61 The boundary conditions 
read 



±i,+i 



(A, C) 



±^poj 2 -^p(t(<s 1 i() -V£(iCICI±iCICI) 



2yfpujiC, T 




Now, we have 



' up 







1 -r <®< 



KKi ± 

1 1 iici 2 

2 ~ 2 Tj 7 



2 



i-C 1 



1±J£ 



with 
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2!) 



and, hence, 



J*M * if a 



, .,3 - 1 - , .3 



icr 



v ^ (u 2 + | C | 3) ._ 4u ,| C | 3 /i«; ! 



1 1 r^tu, if a = — 1. 



This implies 



u: \ y — J— \ ^ 2 + K| 2 1, _ - \ 

i^^e^w — _ /77 A e i, .,2_Li/-i2i2_ A .,i^|3 "to — P^e 



J^ Ae (ur>+K| 2 ) 2 -4w|C| 3 



(^ 2 - ICI 2 ) 

" 2 +ICI 

( w ^-|CP)+4- 



2 +KI 2 



( w 2-|CP)^+4(^-|C| 2 )|CI 2 +4|C| 4 -4^|C| 3 ^ 

|£|2 ^ l(J 



^ 2 +ICI 2 t 

(( w 2_| C |2 )+2 | C |2 ) 2_ 4w | C | a »» 



for a = +1 and 



[p] = k\ s z u 



ICI 1 



ICI 2 



ICI 2 



for a = 
scheme 



T. Therefore, the desired solution may be obtained by solving the splitting 



peu + pd t u — pAu = — Vp, — Ap = in (0, oo) x M™ 
P r B ±1 (u) = 0, [divu] = 0, [p] = S ±1 h w on (0, oo) x 5K™ 
u(Q) = in M" , 
where the bounded linear operators 

5±i : N i(oo) _^ Lp ((0, oo), W^-V^n-i)) 
are defined via their Laplace-Fourier symbols 

rA = ^ 2 ±ICI 2 

i. e. S +1 = (pe + pd t - 2pA T )(pe + pd t - ApSA^ 1 with 

5 := (pe + pd t - pA r ) 1/2 {(pe + pd t - pA T ) 1/2 + (- M A r ) 1/2 } _1 

and 5 _1 = (pe + pd t )(pe + pd t — 2 J uAr)~ 1 , e. g. by the "H^-calculi of the involved 
operators. Indeed, the calculations in Appendix B imply the solution of the above 



splitting scheme to solve the shifted Stokes equations, cf. relation (T^i). On the 



other hand, any solution to the shifted Stokes equations with h w = satisfies the 
equations of the above splitting scheme with a boundary condition 

[p] = hp on (0, oo) x <9R+ 



for some h p € L p ((0, oo), Wp 1 ^ p (W n 1 )). Employing relation (jTi^p once again 



we 



infer h p = S' ±1 (-2^[9 y u;] + [p]) = S* 1 ^ = 0. Therefore, p = and u = 0. 
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5. The Bent Halfspace Case 

In this section, we prove Theorem 12. 31 in the bent halfspace case, i. e. we assume 
a > 0, O = R" with a sufficiently flat function ui G BUC 3 ~ (W^ 1 ) and 1 < p < 
oo with p 7^ |, 3. Following Remark 12.41 it is sufficient to prove existence of a 
unique maximal regular solution for all data satisfying the stated regularity and 
compatibility conditions. However, following Remark l3.7l we may restrict the proof 
to the case a G { -1, 0, +1 }, G { 0, +1 }, / = 0, g = 0, V r h = and u = 0, i. e. 
we have to construct a unique maximal regular solution to the Stokes equations 

pd t u - pAu + Vp = 0, div u = in (0, a) x R^, 

P r B a {u) = 0, QrB f3 {u, p) = h w , on (0, a) x dR™, 

u(0) = in R£, 

where 7i„ := V T h = and h w := Q r /i € o-H"p((0, a), W^ 1/p (r, AT)) PI N^^a), if 
/? = resp. ^ G N^ 1 1/2 _ 1/2j) (a), if £ = +1. 

To solve the above bent halfspace problem, we employ the pull-backs u := tio9 u 
and p := p o 9 W , where 

© w (ar, 2/) := (ar, V + u(x)), x G R"~\ y > 0, 

to reduce the problem to the halfspace case. Since u = u o 0J 1 we have 

d t u = (d t u) o 9J 1 

d k u = (d k u) o 9J 1 - (d k u>) { (dyii) o 9J 1 } 

dlu = (d 2 k u) o e- 1 - 2{d k u) { (d k d y u) o e- 1 } 

+ {d k uf {(dlu) o e- 1 } - (d 2 k u,) { (d y u) o ej 1 } 

dyU = (8yU) O 

dlu = (dlu) o e- 1 

for k = 1, 2, . . . , n — 1 and, hence, 

divu = (divw) o 9- 1 - (Vw| (a„<;) oGj 1 ) 
Am = (An) o 0- 1 - 2( Vw | (V x d y u) o 9J 1 ) 

+ |v^| 2 { (<?^) o e- 1 } - (A W ) { o ej 1 } 

in (0, a) x R™. Here we again decomposed the solution as u — (v, w) and the pull- 
back as u — (v, w) into the first n — 1 components and the remaining component. 
Note, that u will be constructed as a solution to the halfspace problem. Hence v 
constitutes its tangential part and w constitutes its normal part. 

Using this transformation, we first derive the system of interior partial differential 
equations 

pd t v - pAv + V x p = pF v (v, p) in (0, a) xl", 
pd t w — pAw + d y p = pF v (w) in (0, a) x R™, 
divw = G(v) in (0, a) x R£, 
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where; 



pF v (v, p) 



2fj,( Wlu I V x )d y v + n\Wu\ 2 dp - /z(Aw) d v v + (Vw) d. 



pF w (w) := — 2fi(Vuj\V x )dyW + fi\Vuj\ 2 dy-w — {i(Au>) d y w, 
G(v) := (Vw|a w u). 
Now, if we assume HVwH^fj^-i) < S with 6 € (0, 1), we may estimate 

||(Vw| Vx)^,«|| Y/(o) < £|M| Xu(o) , 

HI V -| 2 ^llY /( a)^ 5 ll 5 lk(a), 

ll(Vw)^p|| Y/(o) < <5||p|l x ^i (a) 

and analogous estimates are valid for the terms appearing in the definition of F w . 
Moreover, the embedding chain 

o# p 1/2 ((0, a), L p (R n + , M"- 1 )) ni p ((0, a), H^, K" -1 )) 

^ BUC((0, a), L p (Rl, M"" 1 )) 

^ Y/(a) 

is valid and the corresponding estimates read 

H^lly^a) — al ^ P ll^ S H -B;7C((0,a),Lp(R^,R"- 1 )) 

< ca 1 / p ||a y w|| oH v2 ((0 ^^ Lp(R „ iR „_ 1))nip((0 ^ ) ^i (R „ !R „_i )) , 

where the constant c > is independent of a e (0, 1] thanks to the homogeneous 
initial condition. Hence, 

ll(Aa;)^t;|| Y/(o) <ca 1 /i'||Aa;|| i;oo(RB _ 1) ||t;|| Xti(o) 

and an analogous estimate is valid for the term appearing in the definition of F w . 
Summing up the above estimates, we derive 

(29a) * )IIy ' (o) " M3<5 + caVP ll Aw llioo(R«-^)H"llx tl (a) + S \\P\\xi^ a y 

ll^(«;)|| ¥/(o) < M3 ( 5 + ca 1 /p||A W || L3o(R „_ 1) )||z7)|| Xii(a) . 

Furthermore, we have 

ll( Vw l^)llo^ a ((0 1 a),L p (R ¥ ))^ 5 ll^llo^ a ((0,a),i p (R ? )) 

and 

\\(Vu\d v v)\\l MOa)Hlp{Rn+)) = \\(^\dyv)\\ P Lp((aa)LpiRV) 



n-l 

I a. V7, , I £> -, Ml 

Lp((0,o),L p (R!f.)) 



fe=l 
n-l 



J2\\(^\d k dyV] 



lll,p((0,a),Lp(M™)) 
fe=l 



+ ll(V-|^)t ((0 , a) , ip(R n )) 
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with the same arguments as above implies 

|| ( Vw | d y v )|| Lp((0 , a) , W)) < 2 1/P (S + caV'W VMl^-^Mx^ 
where the constant c > is again independent of a £ (0, 1]. Hence, the estimate 

(29b) \\G(v)\\Y g (a) < 2 1/P (<5 + ca^\\VM L ^-i))\\ 

is valid. 

Now, the outer unit normal of R™ is given as 



vr(x, uj(x)) = 



1 



lVw| 



Vw(ar) 
-1 



X £ 



*(<*) 



pn— 1 



and we infer 



{ VY } = ^£ + 



Vw 




=: Vy, + 7V(Vw), 



where we have set k w := y 1 + || Vw|| and v^, denotes the outer unit normal of K™ 
on E := Moreover, 

' Vw(x) ® Vw(x) -Vw(i)' 
-Vw(x) T 1 



Pr(x, u(x)) = 1 5- 

1 + llv^H 2 



X £ 



jn- 1 



and we infer 



Pr.- 



1 + |Vw| 



Vw ® Vw -Vw 



P s - L(Vw) 



-Vw -|Vw| 

If a = 0, the tangential boundary condition Pr[tt] = on (0, a) x <9R™ is equiv- 
alent to 

Pe[U] - L(Wuj)[u] = P r [u] o e w = 
on (0, a) x QR™, which may be rewritten as 

[v] = L„ ( Vw) [u] , L„ (Vw) [«] = 

on (0, a) x cM™ , where L„ denotes the first n — 1 rows of L and L w denotes the 
last row of L. However, the first equation already implies the second one, which is 
not surprising, since a tangential boundary condition on T = dQ may not lead to n 
linearly independent transformed boundary conditions on E = 9R™. The original 
boundary condition in the bent halfspace is therefore equivalent to 

[v] = H%(v, to) on (0, a) x OR^ 

where 



H° v {v, w) :-- 



{(Vw\[v])-[w]}. 



l + ||Vw|| 

Observe, that all terms above either carry a factor Vw or are of lower order. With 
the same arguments as above we therefore infer 

(29c) \\H°(v, w)\\ K(a) < c^ + o^HV^II^-^)) 11 --- 11 

with some constant c > independent of a £ (0, 1]. 



*(«) 
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If a = ±1, we have to treat the terms 
Vu ± Vu T ] o9,= Vw ± Vu T ] - 

[ ±(Vw|[fyt>]) 

=: [Vu ± Vu T ] - M(Vw, [fyv], 

and the tangential boundary condition ^Pr[Vw ± Vu T ]^r = on (0, a) x <9R™ is 
equivalent to 

uP^Vu ± VS> S - uL(Vlu) { [Vu ± Vti T ] oG,}^ 

- /xP E M(Vw, [fyt;], [a y u)])i/ s 

+ /i{Pr[Vu±Vu T ]o9 w }iV(Va;) 
= k^M { Pr[Vu ± VuV o 6 W } 

= o, 

on (0, a) x <9R™ , which may be rewritten as 

T£t[dj/v] _ MVxH = fiL v (Vu}) { [Vu ± Vti T ] o6 u }i/ s 
+ uM v (Vu, [d y v], [d y w])vy, 

- u{P^[Vu± Vm t ] o 9 W } 7V(Vu;) 

and 

uL w (Vlo) { [Vu± Vw T ] o9 u }^ = ^{ Pr [Vm± Vw t ] o 9 w }iV(Vcj) 

on (0, a) x <9R™ , where M„ denotes the first n — 1 rows of PeM, Pp denotes the 
first n—1 rows of Pr and Pp" denotes the last row of Pp. Again, the first equation 
already implies the second one and the original boundary condition in the bent 
halfspace is therefore equivalent to 

Tuidyv] - uV x [w] = H^(v, w) on (0, a) x dW* 

where 

H±\v, w) := uL v (Vuj) { [Vu ± Vu T ] o0 u }^| uM v (Vu, [d y v], [d y w\)vr. 

- n { P£ [Vu ± Vu T ] oQ u }N(Vlu). 

Observe, that all terms above again cither carry a factor Vu or are of lower order. 
With the same arguments as above we therefore infer 

(29d) \\H^(v, w)\\ TW < c(<5 + a 1/p ||V 2 o;|| sc/cl _ (R „_ 1) )||tZ|| Xit(a) 

with some constant c > independent of a e (0, 1]. 

Concerning the normal boundary condition, we start with the case (3 = 0, i. e. 
[u] ■ vr = h ■ vr on (0, a) x <9R™. This condition is equivalent to 

-[w] + Vlo ■ [v] = [u] ■ vy, + [u] ■ N(Vlj) — { ([u] ■ v T ) o 6 W } = K u h w 
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on (0, a) x <9R™ with h w := (h-ir)oQ u . Therefore, the original boundary condition 
in the bent halfspace is equivalent to 

[w] = Vw • [v] - K u h w =: H°(v) - n u h w 

on (0, a) x As above, all appearing terms either carry a factor Vu; or are of 

lower order and we infer 

(29e) \\H° w (v)\\ K{a) <c(J + a 1 /P||V 2 w|| Bt , cl _ (RB _ 1) )||ti|| Xu(o) 

with some constant c > independent of a e (0, 1]. 

Finally, if j3 = +1, the normal boundary condition 2[id v u ■ vr — \p] = h ■ z/r on 
(0, a) x cM™ is equivalent to 

2fj,d v u ■ vt. - [P] - 2[j,Q(Vlo, [d y v], [d v w])v?, ■ i>t. 

+ 2/z {[Vu] o 9 w }iV(Vu;) • z/ E 
+ 2[i {[Vu]v r o 6 W } • K w iV(Vu;) 

- \^\ 2 M 

= K l {{2iid v u ■ vt) o9 u - \p] o O^} 

on (0, a) x <9R™ with h w := h- fr 9 W and 

Vw <g> (Vu) [<9 y w] 



[Vu] o6 u = [Vu] 



=: [VS]-Q(Vw, [a y tJ], [0„ti>]). 





Therefore, the original boundary condition in the bent halfspace is equivalent to 
2[id v u ■ v-z - [p] = 2fiQ(Vuj, [d y v], [d y w])u^ ■ vy, 

- 2n{[Vu}oe LU }N(Vu J )-vz 

- 2 M {[Vu]^ r oe4-K U 7V(Vw) 

+ \Vuj\ 2 [p] + n 2 Ji w 

= ■ H+ X {v, w, p) + n 2 Ji w 

on (0, a) x dW± ■ As above, all appearing terms either carry a factor Vw or are of 
lower order and we infer 

(29f) \\H+\v, w, p)\\ K{a) < c(S + o 1 /f||V 2 W || Btrcl _ (RB _ 1) )(||«|| Xu(o) + \\P\\ <y{a) ) 

with some constant c > independent of a € (0, 1]. 

Summarizing the above considerations, the desired maximal regularity property 
of the bent halfspace problem is equivalent to the existence of a unique solution 
(v, w) e o#p((0, a), L p (Wj_, R n )) n L p ((0, a), , W 1 )) together with an ap- 

propriate pressure p, whose regularity depends on j3 and 7, to the fixed point 
problem 

L a (v, to, p) = R a (v, w, p), 
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where L a resp. R a denote the left resp. right hand side of the transformed system, 
i. e. L a is the isomorphism arising from the halfspace problem. Now, the estimates 
([29)) imply, that we may first choose S E (0, 1) and then, given w g BUC 3 ~ (R rl_1 ) 
with |Vw| L (R"-i) < $i we may choose a € (0, 1], such that the fixed point problem 
admits exactly one solution (v, w, p) in the maximal regularity class, thanks to the 
contraction mapping principle. Hence, the original problem in the bent halfspace 
has the desired maximal regularity property at least on small time intervals. Ob- 
serve, that a may depend on uj, but not on the data. Therefore, we may construct 
a unique solution on any arbitrary time interval by successively solving the system 
on sufficiently small time intervals. Finally, the bent halfspace problem has the 
desired maximal regularity property for all a > and all w € _B[/C 3 ~(R™~ 1 ) with 

6. The Case of a Bounded Domain 

In this section, we prove Theorem l2.3l in the bounded domain case, i. e. we assume 
a > 0, ft C M" to be a bounded domain with boundary T — dtt of class C 3 ~ and 
l<p<oo,p^|,3. Following Remark 12.41 it is sufficient to prove existence of 
a unique maximal regular solution for all data satisfying the stated regularity and 
compatibility conditions. However, following Remark l3.7l we may restrict the proof 
to the case a £ { —1, 0, 1 }, j3 € { 0, 1 }, / = 0, g = 0, Prh = and u = 0, i. e. we 
have to construct a unique maximal regular solution to the Stokes equations 

pdtu — fiAu + Vp = 0, divu = in (0, a) x f2, 

(30) P T B a {u) = 0, QrB f3 (u, p) = h w , on (0, a) x dQ, 

u(0) = in Q, 

where h v := P v h = and h w := Q v h e N^ 7 with h w € o#p((0, a), Wp 1/p {T, NT)) 
and 7 = -co, if /3 = resp. 7 = 1/2- l/2p, if /3 = 1. 
If f3 = 0, we first choose an offset 

V77 e H^((0, a), L P (Q, R n )) n L p ((0, a), i/ p 2 (fi, E n )) 

as a solution to the elliptic problem 

— Ar] = in (0, a) x Q, d v r\ = h w ■ v on (0, a) x <9f2, 

we set h v := — PrS Q (V?y), and then construct a unique maximal regular solution 
(u, p) := (u — p + pdtrj — liArj) to the Stokes equations 

pdtu — /.tAu + Vp = 0, divu = in (0, a) x il, 

(31) P r B a (u) = h v , QtB^u, p) = 0, on (0, a) x Oil, 

u(0) = in Q, 
where we dropped the bars again. 
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Now, if we consider (f3lj) in case ft = and (f30|) in case j3 = 1, any maximal 
regular solution, where we may assume 



(32) 



if ft = 0, enjoys the additional time regularity 

p G „££((Q, a), ^p(fi)) n L P ((0, a), ff^Q)) 

for all (7 G (0, 1/2 — l/2p), which can be seen as follows. Given ip G where 
1/p + l/p' = 1, we choose G H 2 ,{VL) to be a solution to 

— = in ^> = on <9S7, 

if /3 = 0, where ipo := ip — resp. 

-A(f> = ip in Q, [<j)} = on dfl, 

if /3 = 1. Observe, that (|3"2"j) implies ( p ] V'o ) = (P I V0> if /3 = 0- Using integration 
by parts and the interior equations we infer 



(p\i>) = -(p|A0) = (Vp|V^)- J[ P ]d^da 

r 

= -pd t (u\V<l>) + iJL(Au\V(l>)- J\p]d u <l>da 

r 

= /i( Au| V0) - p9 f / ([it] • i/) [</>] do- - / [p] a„0dcr 



and the boundary conditions imply 

(p\ip) = n(Au\V<t>), 

if /3 = 0, resp. 

(p|V) = A»( A «|V0)- / gd^dcr, 



if /3 = 1, with 5 := ( h w \ v ) — 1[i d v u ■ is, i. e. 

q G % 1/2 - 1/2p ((0, a), Lp(50)) n L p ((0, a), ^"^(flO)). 
Employing an integration by parts again, we have 

H( Au | ) = M J d v u ■ [V0] dcr - n( Vu \ V 2 (j) ) 



and we infer 



(p\tp) = fj, J d v u- [V0] da - J qd u (f)&a - Li{Vu\V 2 (j)), 
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where we have set q := 0, if ft — 0. Now, we may use the regularity of the functions 
involved on the right hand side and apply the operator d% to obtain the estimate 

H^tPll.Lp((0,a)xfi) 

< c (||9fa^|| Lp((aa)xr) +||af g || Lp((0a)xr) +||afVu|| Lp((0a)xn) ) 

and, hence, the desired regularity property of the pressure p. To be precise, we 
have the estimate 

lblU ? ((0,a),i p (O)) < C (lMU„(a) + PW h W Ytf(a)) ■ 

Note, that for an inhomogeneous right hand side / € oHp((0, a), L p (ft)) in the 
momentum equation of (|30j) respectively (|31|) . an analogous computation yields 

\\p\\ H°{(0,a),L p (n)) < C (lMI X u (a) + \\f\\ H°((0,a),L p (n)) + ^INU^(a)) ' 

Note, that the embeddings 

i?;/ 2 ((0, a), X), oH;((0, a), X) L g ((0, a), X) 

are available for some p < q < oo and every Banach space X, where the embedding 
constant does not depend on a > thanks to the homogeneous initial condition. 
Hence, we always have 

i?; /2 ((0, a), X), q H;((0, a), X) L p ((0, a), X), 

with embedding constant ca r , c > being independent of a > 0; r = 1/p — l/q. 

To prove the existence of a unique maximal regular solution to the Stokes equa- 
tions pop resp. (|3ip we first choose finitely many points x\, x%, ■ • • , xn G T, such 
that r n B r (xk) is the graph of a BJ7C 3_ -function cjfc over the tangent plane T Xk T 
for fc = 1, 2, JV. Moreover, we choose r > sufficiently small, such that 
||Vwfc|| L (Rn-i) < S for k = 1, 2, . . . , N with S € (0, 1) as in the previous section. 
The open sets ?/& := B r (xk) then constitute a covering of T, which may be com- 
pleted by an open set Uq C to a covering of f2. Finally, we choose a partition of 
unity 4>o, fa, . . . , </>jv € Co°(R n ) subordinate to the covering Uq, U\, . . . , Un of fl 
and cut-off functions V'o, ipi> ■ ■ ■ , "4>n G C^°(IR n ) with spt^fc C [/& and t/j^ = 1 on 
spt <?!>fc. 

Now, every maximal regular solution (u, p) of system (|30[) rcsp. (|3ip may be 
decomposed as u = no + u\ + ■ ■ ■ + un and p = pa + pi + ■ ■ ■ + pn with Uk = 4>kU 
and pk = 4>kP for k = 0, 1, . . . , N. By construction, (uq, po) is a maximal regular 
solution to the whole space Stokes equations 



pd t u a — /uAu + Vpo = pFo{u, p), divM = Gq(u) in (0, a) x R", 

u (0) = inE n , 
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and (uk, Pk) is a maximal regular solution to the bent halfspace Stokes equations 
pd t u k - pAu k + Vp k = pF k (u, p), d\vu k = G k {u) in (0, a) x 
P rk B a (u k ) = H« k (u) + h Ttk , 
Qr k B (u k , Pk) = H^ k (u) + h v , k 

Ufc (0)=0 inR£ fc 

for k = 1, . . . , N with 



on (0, a) x 



pF k (u, p) 
G k (u) 

H T,k( U ) 

KM 



= -n[A, <f> k }u + [V, 4> k ]p = -p(A<f) k )u - 2/i(V<fo • V)« + (V0 fc )p, 
= [div, = V^fc • u, 
= [Pr fc ^°, - 0, 

= [Pr^ 1 , <p k ]u = pPr k (V0 fe <g> u ± u <g> V(p k )v k , 

= [Qr k B°,^ k }(u) = 0, 

= [Qv k B l , (j>k\{u) = 2p([V(j) k ® u]v k ■ Vk)v k 



and h T _ k = 4> k h v , h v _ k = 0, if j3 = 0, resp. h T) k = 0, = <j) k h w , if /3 = 1. 
Therefore, w = v Q + V770, Po = <?o - P<9 t ?7o + ^Ar/o, and u k = v k + V-q k + v k , 
Pk = qk~ pdtVk + P^Vk + q~k for k = 1, . . . , N, where 

770 G o^((0, a), H 2 p (R n )) n # p 1/2 ((0, a), ff p 3 (R")) 

is a solution to the wholespace problem 

-A770 = -G Q (u) in (0, a) x M™ 

and 

Vk € H*((0, a), H*(R2 k )) n H^dO, a), ff p 3 (R£J) 
is a solution to the bent halfspace problem 

-A Vk = -G k (u) in (0, a) x R" fc 
complemented by the boundary condition 

d v t] k =0 on(0, a)xdSZ k , 

if /3 = 0, resp. 

[»fc]=0 on(0, a)xdSZ h , 

if /3 = 1 for k = 1, . . . , N. Moreover (vo, Po) is a maximal regular solution to the 
wholespace Stokes equations 

pd t v - pAv + Vgo = pFo(u, p), divu = in (0, a) x R™, 
w (0) = inM™, 
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{vk, Qk) is a maximal regular solution to the bent halfspace Stokes equations 

pd t v k - pAv k +\7q k = pF k (u, p), divu fc =0 in (0, a) x R" fc , 
Pv k B a {v k ) = H« k (u) - P rk B a (V Vk ), 
(34) Qr k BP(v k , q k ) = H^ k {u) on (0, a) x cM» fc , 

- Qr k B p {\7r] k , pAr] k - pd t r) k ) 

v k (0) = inMZ k 

and (v k , q k ) is a maximal regular solution to the bent halfspace Stokes equations 

pd t v k - pAv k + Vq k = 0, div v k = in (0, a) x , 
Pr k B a (v k ) = h T , k , 
Qr k B^(v k , q k ) = h v<k 

v k (0) = in 



(35) _ on (0, a) < . 



for k = 1, . . . , N. Now, every term on the right hand side of (|33|) and (134)) carries 
additional time regularity. To be precise, we may exploit the extra time regularity 
of the pressure to obtain 

\\pF k (u, P)\\ oH? ((0,a),L p (M.n k ,K»)) ^ c (|I m II X„(q) + \\p\\ H^ ((0,a),L p (U))J 

^ C (lMU„(a)+/?ll^ll Y^(a)) 

and, therefore, 

\\pF k ( U , P)|| Lp((0 ,a)xR Sfc ,K") < CaT (ll U H X„(a) +/ 3 ll /l ll Y^(a)) 

for k = 0, . . . , N . The commutators of the boundary conditions are of lower order 
and we have 

\\Kk^) + KM\\ aY ^ {a) <^ T \\u\\ M 

for k = 1, . . . , N, If P = 0, we have Qr k B l3 (Vr] k , pArj k — pd t rj k ) — and, thus, 

IMI X u (a) + IMIx^a) ^ Ca l U HoX u (a) 

as well as 

Wlk\\ oH? ((0.a),L p (RZ k )) ~ ° {^ Vk hxu(a) + \\P F k( u > P) \\ o H ? ((0,a),L p (M» k , R»)) 
^ C IMU„(a)- 

for k = 0, . . . , N. If p = 1, we have 

Qr k B p (Wr] k , pAr] k - pd t r\ k ) = 2p\^ 2 r\ k \v ■ v - [pAij k - pd t r] k } 

= 2p[S7 2 Ti k }v ■ v- p\Ar] k \. 
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However, 

H V2?? ' i; lloHp /2 ((0,a),L p (R Sfc ))nL p ((0,a),Hi(R Sfc )) 

< ca r ||V 2 w fe || o ^ ((0 Q) Lp(R „ fc))nff i/2 ((0 a) iff i (R „ fc)) 

< CaT H U IU„(a) 

and we again infer 

IMU„(a) + ll^llx^fa) < CaT |l M HoX u (a) 

as well as 

hk\\ a H ? ((0,a),L p (M2 k )) ^ C (KlloX,(«) + IIP**(«> V) \\ Q H° ((0,a),L p (R» fc , R»))) 

+ c||[V 2 % ]|| oY ^ (a) 

^ C IMU„(a)- 

for k = 1, . . . , N. Finally, we have 

INU„(a) + ll%llx^(a) < c||/l|| oY ^ (a) 

and 

llftlll W a),L p (R» b )) < C (lNU(a) + ^INU^a)) < C INU^(a) 

for fc = 1, . . . , AT. Now, these estimates imply 

I^Aryo - pdtVo\\ oH; ao.a).L p (R Zk )) = \\Po - Qo\\ oH? ((o,a),L p (mz k )) 

<c(NU(.) + W ^ w ) 

and 

H^Affc - pft»?fc|| oH?(( o,„) 1 i, J( (R» (i )) = HP* - 9fc - 9fc|| ff ? ((0 1 o),£p(R» Ji )) 

^ C ( l|u|l o^W + l|/l|l o^(a)) 

for k = 1, . . . , N and since 

ll A ^IU p ((0,a),L p (R Sfc) ) < C IMI X„(a) 

we further infer 

\\pdtVk\\ oH?ma ),L p (m Zk ))' \\^Vk\\ oH?m a),L p (MZ k )) ^ C (lMU„(a) + H^U^a)) 

for k = 0, . . . , N. Thus, using 

ll^ ? 7fe|lL p ((0,Q),ff|(E S)c )): ll A? ?fc|| Hi((o,a),L !) (K 3fe ))nI, p ((0,o),H|(K 3fe )) < C ll U ll X„(o) 

we also have 

W d ^Vk\\ LpmahLp{M , k)) < caT \\ d tm\\ QH ^ mahHl( ^ k)) 

< ca T \\dtVk\\ oH?ma) ,i p (R» ji ))ni, p ((o,o),ff»(H» lb )) 

< ^IblloX^^+cll^lloY^Ca) 
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and 

W^Vk\\ LpmahLp{MZk)) < «* T H A »Zfc|l Hi/ a ((0,a) > Hi(R» fc )) 

< ca T ||A77 & || off i ((0ia)iip(K ^ ))nip((0)0)iH | (K „ fe)) 

^ CaT H U H X„(a) 

for fc = 0, . . . , N. Hence, we have 

ll V7 ? fe ll X„(a) + \\^ A Vk ~ Pdtm\\^ n(a) < CaT H M ll„X„(a) + C INI Y^(a)' 

which together with the above estimates implies 

IKH X u (a) + bfellx^(a) ^ CaT |MU„(a) + C INI Y^(a)> 

Finally, we have 

N N 

(«> p) = ^2(^ku k , ipkPk) = ^2(u k ,p k ) 

k=0 k=Q 

and choosing a > sufficiently small, we deduce the uniqueness of the solution as 
well as the maximal regularity estimate. 

To prove solvability of the Stokes equations, let a € { —1, 0, 1 }, j3 € { 0, 1} and 
7 = — oo, if f3 = 0, resp. 7 — 1/2 — l/2p, if /3 = 1. We consider the bounded linear 
operator 

L Q,/3 : o^(a) — ► Y"'^(a) 
defined by the left hand side of the Stokes equations, where we set 
X^(a) := { (u, p) G X u (a) x x£ 7 (a) : u(0) - } 

and define 

Y«^(a) := { (/, ff , ft) € Y f (a) x Y» x Y^(a) : (/, <?, ft, 0) € Y«^(a) } . 

By the above considerations, L Q -^ is injective with closed range and it remains 
to prove surjectivity. To accomplish this, it is convenient to construct a bounded 
linear right inverse 

: Y^(a) — > X?(a). 

First, let j3 = and 7 = —00. By Theorem 13.61 and the construction at the 
beginning of this section, there exists a bounded linear operator 

S?' : oY^(a) — ► „X*(a), 
such that (w, p) = S"' (f, g, ft) satisfies 

pdtv — fiAv + Vg = pf, div v = g in (0, a) x f2, 
Q T B p (v) = Q r ft, on (0, a) x dQ, 
v(0) = in a 
Now, we construct a bounded linear operator 

S?" ■■ rY£>) := { t, € Y$(a) : Q r?7 = 0, r?(0) = } — > X» 
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and then set 

S a >?(f, g, h) := S?' (f, g, h) + S^(P r h - P r B a S^(f, g, h)) 

for (/, g, h) e Y"'P(a). Given h e T Y^(a), we define (v k , q k ) to be the solution 
to the bent halfspace Stokes equations 

pd t v k - pAv k + Vq k = 0, divu fe =0 in (0, a) x M" fc , 

Pr k B a (v k ) = <t> k P T h, Q Tk BP(v k ,q k ) = on (0, a) x dS£ k , 

v k (0) = inM2 fc 

for k = 1, . . . , N. Then we define r) k G #p((0, a), if* (ft)) n Hp /2 {{0, a), 
to be the unique solution to 

— Ar] k = —Vip k ■ v k in O, d v r\ k = on 90 

again for k = 1, . . . , N. Finally, we set 

N 

{v, q) = ^OfcW/c - V%, q k + pd t T] k - p,Atj k ) 
fe=i 

and define S^'^h :~ (v, q). This way for (/, g, h) e Yy^(a) we have 
L a,PS°>P(f, g, h) 

= L^S^(f, g, h) + L^S^(P r h - P r B a S^(f, g, h)) 
= L a 'P(v, q)+L a 'P(v, q) 

with 

(v, q) = S^(f, g, h), (v, q) = S^(P r h - P r B a v). 

Now observe, 

N 

pd t (v + v)- /j,A(v + v) + V(q + q) = pf + Y,(- f*l A > V>fcH + [V, Vfc]%) 

fe=i 

by definition of and S% Ji . Furthermore 

N 

div (v + v)= g, P r B a (v + v) = P r h + Y / ([Pr k B a , VfcH - P r B a (V Vk )) 

fc=i 

as well as 

Q T B p (v + v, q + q) = Q T h and (v + v)(0)=0. 
Hence, L a ^S a ^ = I + R a - p with 

N 

Rf {f, 9, h) = ^2 ( - V'fcH + [V, Vfe]%) 
fe=i 

and g, /i) = 0, QrR^if, 9, h) = as well as 

AT 

Pr^'^(/, 5, M = E ([ p r fe ^ Q , V>fcH - W(V%)). 
k=i 
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The subscripts denote the different components of R a 'P . Now, 

IMU u (a) + Wlk\\ X ^ {a) < C ll^ll Y^(a) 

with h = P r h- P r B a v and 

hk\\ oH? ((0,a),L p (^ k )) ~ C IHL„(a) < C IWI Y^(a)' 

which implies 

|| [V, Vfe]ft-A*[A, 1pk]v k \\ Yf (a) < CflT (lHI X„(a) + lkfc|l o ^((0,a),L p (M 3fe ))) 

< ca T \\h\\ oYl:>) , 
since the commutator terms are of lower order. Analogously 

\\[Pr k B a , ^kW\\ Q Y^(a) ^ CaT H W fcll X u (a) < ^ 1 1 M 1 Y^ (a) ■ 

To estimate the final term PrB a (Vr] k ) we set v k ■= 4>kVk — Vr/^, q k ■= tpklk + 
pdtilk — IJ>Ar]k and observe, that 

pd t v k - pAv k + Vq k = -/z[A, ip k ]v k + [V, ^k]qk, divw fc =0 in (0, a) x fi, 
P T B a (v k ) = <p k P r h + [P Tk B a , VfcH - P T B a (Vr)k), 

on (0, a) x 90, 

Qr k BP(v k , g k ) = 

w fc (0) = in O. 

Now, 

||Wfc|| oXu (o) ^ C (\\ V k\\ X u (a) + ll9fc|li p ((0,o)xR» fc )) + ll V %H X„(a) + H^H Y^(o)) 

and 

|| [V, lP k ]q k - fl[A, *l>k]Vk\\ H*UO,a),L„pZ k )) 

' (\\ V k\\ X u (a) + Wlk\\ H ? ((0,a),L p (R Zk)) ) < c||ft|| o¥ -.fl (o) 



imply 



hk\\ oH? ((0,a),L p (RZ k )) 

+ ||[V, tpkhk - /i[A, ^fe]Wfe|| H-((0,o),L p (M 3fc ))) 

< c||ft|| oY -.fl (o) 

and we infer 

||/iA»fc - M%|| oi ^ ((0 , a )„L p (fi)) = HVfeSfe - 9fc|| 0ff? ((0, o ),L p (n)) ^ C IWI Y^(a)- 

Since 

H A %IU p ((0,a),L p (Q)) < C KH X u (a) < C ll^ll Y^(a) 

we further infer 

\\pdtVk\\ oH?mahLpm , \\^Vk\\ oH?ma) , Lpm < c||ft|| oY a,p (o) . 
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Thus, using 

ll a «wl£ p ((0,a),H2(n))> \\ A Vk\\ oH i ({0 ,a),L p (n))nL p ((0,a)M^(n)) 

<c||Wfe|| oXu(a) < C||^|| o¥ a,j» (o) 

we also have 

\\dtVVk\\ Lp ((0,a),L p (n)) < mT \\ d M\ oH °/* {{0ia)tHim 

< ca T \\d t r]k\\ oH a {(0 ^^ Lp(n))nLMO a) m {Qj)) 



< ca T \\h\\ 



and 

W^Vk\\ LpmalLpm < ca-||A % || oH x /2((0i0)i ^ (n)) 

< ca T || ^Vk\\ oH i^o,a),L p (Q))nL p ((0,a),H}(Q)) 

< ca-fh\\ aK: , {ay 

Hence, we have 

\\Vrik\\ oXu{a) + W^Vk - p9tVk\\ x nia) < ca T \\h\\ aY ^ (a) , 

which implies 

||P r B«(Vr ?fc )|| o¥ ^<c^||ft|| o¥ ^ (o) . 

Since 

IWU^(a) ^ c ll(/> 9, h)\\ aY ^ (a) 

we obtain 

\\R a ' P {f, 9, h)\\ BY ^ {a) <ca^\\(f,g, h)\\ oYr(ay 

Thus, choosing a > sufficiently small, R a ^ is invertible by a Neumann series, 
which yields the right inverse S a ^(I + R 01 ^)' 1 for L a ^ . 

The case (3 — 1 and 7= 1/2 — l/2p may be treated analogously. The only 
difference is to use a Dirichlet problem to construct the r)k instead of a Neumann 
problem. However, a more elegant proof is also available by a homotopy argument. 
In fact, one may repeat the arguments of Sections @] to [HI to infer, that the bounded 
linear operators 

defined by the left hand side of the Stokes equations 

pdtU — /iAu + Vp = pf, divu = g in (0, a) x fl, 
PrB a u = Prh, t 2p d v u ■ v — [p] = h ■ v, on (0, a) x <9f2, 
u(0) = in n. 

are injective with closed range for all t 6 [0, 1]. Hence, these operators are semi- 
Fredholm. Thus, Lf = L a - +1 is Frcdholm of index zero and, hence, an isomorphism, 
since Lq = enjoys this property by Theorem 13.61 This completes the proof 

of Theorem [ 
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7. Local Well-Posedness of the Navier-Stokes Equations 

Finally, we prove Theorem 12.11 for a bounded domain SI C M. n with boundary 
r := dVl of class C 3 ~. Following Remark 12. 2\ it is sufficient to construct unique 
local-in-time strong solutions to the Navier-Stokes equations 

pd t u + div(pu eg u — S) = pf, divu = in (0, oo) x Q, 

(N)^' B B(u, p) = on (0, oo) x dQ, 

u(0) — uq in Q 

where B = B a 'P with a, f3 E {— 1,0, 1} realizes one of the boundary conditions 
(B), for given data / G £ p ((0, oo) x f2) and uq G VFp 2 ^ p (£l), where we assume 
rt + 2 < p < oo. 

First, let (it*, p*) be the unique maximal regular solution to the Stokes equations 
(S)/o'o«o' Then, P) := ( u — u * ' P — P*) ^ s a solution to 

pd t u — fiAu + Vp = N*(u), divtt = in (0, a) x ft, 

B(u, p) = on (0, a) x dfl, 

u(0) = in ft 

for all a £ (0, 1], where the nonlinearity on the right-hand side is given as 

N*(u) = -((ft + u*) ■ V)(u + u*). 

This is the well-known nonlinear perturbation, which always occurs, if the Navier- 
Stokes equations are reduced to the Stokes equations. Since 

u G o^((0, a), L p (ft, R")) n L p ((0, a), ff p 2 (ft, R")) 

and p > n + 2, we may employ the standard estimates of the nonlinearity N* to 
choose a G (0, 1] sufficiently small, such that existence and uniqueness of a local- 
in-time solution follows by a contraction mapping argument. We do not want to 
repeat these well-known arguments here and close the proof of Theorem 12.11 

Appendix A 

Parabolic Systems with Prescribed Divergence: L p -Maximal 

Regularity 

In this appendix, we establish the maximal regularity property of the parabolic 
problem 

pdtu — pAu = pf, in (0, a) x ft, 
(PT/luo = h on (0, a) x 9ft, 

u(0) = uq in ft 

with a G {—1, 0, 1}, where the boundary condition is given by the linear operator 
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resp. 

Btl( u ) ■= M-Pr(Vu ± Vu> + [divu]u 

for a = ±1. We assume a > and 11 C R" to be a halfspace, a bent halfspace 
or a bounded domain with boundary T := dfl of class C 3 ~ and 1 < p < oo with 
p 7^ |, 3. The claimed maximal regularity property of the system (P)^ is 
therefore equivalent to the existence of a unique maximal regular solution 

u G f£((0, a), L p {ft, R n )) n L p ((0, a), P 2 (i2, R n )), 

whenever 

/ G L p ((0, a), L p (fi, R n )), h G L p ((0, a), L p (r, R™)) and u € Wp 2 " 2 / p (Q, R n ) 
with 

P r h G Wp- 1/2p ((0, a), L p (r, IT)) n L„((0, a), W^ 1/p (T, IT)), if a = 0, 
Pr^ G W p 1/2_1/2p ((0, a), L P (T, IT)) n L p ((0, a), W^ 1/p (T, TT)), if a = ±1 
and 

Qrh G Wp 1 / 2 - 1 /2p ((0 , a), Lp(r, AT)) n L p ((0, a), W^(T, NT)) 
satisfy the compatibility conditions 

P r [u ] = P r h(0), if a = and p > |, 
2pP r (Vu a ± VuJ)z/ = Pr/i(0), if a = ±1 and p > 3 

and 

[divuo] = ft(0) • v, if p > 3. 
We will prove the claimed maximal regularity property of the system (P)^ h uo in 
the following subsections for the different cases a = and a = ±1. 

A.l. The Case a — 0. We treat the system like the Stokes equations in the pre- 
vious subsections by a localization procedure and start with the halfspace case 
= R". We decompose the desired solution as u — (v, w) into a tangential part 
v : (0, a) x M™ — ► R" -1 and a normal part w : (0, a) x M™ — ► R. The system 

(P)-; 

h,uo then reads 

pd t v - fj,Av = pf v , pd t w - pAw = pf w , in (0, a) x M™, 

[v]=/i„, div x [w] + [<9 y w] = on (0, a) x <9R^, 

u(0)=u o , w(0) = w o inR™, 

where we decomposed / = (/„, f w ), h = (h v , h w ) and uq = (vo, u>o) into a tangen- 
tial part and a normal part analogously to the solution u. Moreover, we decomposed 
the spatial variable into a tangential part x € R" _1 and a normal part y > 0. Of 
course, div-r denotes the divergence w. r. t. x. Obviously, this system may first be 
solved for a unique maximal regular solution v, which may then be used as part of 
the data of the parabolic boundary value problem for w. Hence, w may in a second 
step be obtained as the unique maximal regular solution to this remaining problem. 
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Now, we assume ft = R™ to be a bent halfspace with a sufficiently flat function 
ui G BUC Z ~ (R™ -1 ). We may first solve the plain Dirichlet problem 

pdtu — p,Au = pf, in (0, a) x f2, 

[u] = P v h+ (e' Ar [u ] • on (0, a) x 90, 

«(Q) = uq in f2, 

which forms a system of n decoupled parabolic initial boundary value problems 
with Dirichlet condition. Here, Ap denotes the Laplace-Beltrami operator on T. 
Hence, we may in the sequel assume / = 0, Prh = and uo = 0. Therefore, we 
may repeat exactly the same arguments as in Section [5] to reduce the problem to 
the halfspace case. 

Finally, we assume f2 C R™ to be a bounded domain with compact boundary 
r = dQ of class C 3 ~. Since we may solve the above plain Dirichlet problem, we 
may again assume / = 0, Prh = and Uq — 0. Therefore, we may repeat exactly 
the same arguments as in Section [6] to reduce the problem to the whole space case 
resp. the bent halfspace case. The treatment of the case a = is therefore complete. 

A. 2. The Case ex — ±1. Here all boundary conditions are of order one. There- 
fore, the system fits into the context of the abstract parabolic problems, which are 
completely treated in [TU] . Indeed, the system satisfies all necessary and sufficient 
conditions, which reveal L p -maximal regularity - especially the Lopatinskii- Shapiro 
condition. For more details we refer to [SJ [TU] ■ 

Appendix B 

Parabolic Systems with Prescribed Divergence: Special Solutions 
In this appendix, we study the parabolic system 

peu + pd t u — pAu — —Vp in (0, oo) x R", 
(P)-vio Q o B&M = on (0, oo) x dWl, 

u(0) = in R™ 

with a G { — 1, 0, I}, where the boundary condition is defined as in the previous 
appendix, i. e. it is given by the linear operator 

^div(w) ~Pr[u] + [diYu]u 

resp. 

£±i(V) : = f j,P r (Vu± Vu T )v + [divw]^ 
for a = ±1. Moreover, the parameter e > denotes an arbitrary shift. We will 
limit the following considerations to the special case, where the right hand side 
— Vp G L p ((0, to), X P (R™ , R n )) is obtained by either solving the Dirichlet problem 

(E D )^°' R + -Ap = in (0, to) x R™ , [p] = h on (0, to) x dW'l 

for a given function h G L p ((0, to), Wp 1 ^ P (9R")), or the Neumann problem 

(Eat)^ R + -Ap = in (0, to) x Rl, d v p = h on (0, oo) X 
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for a given function h G L p ((0, oo), Wp 1/p (<9R™)). Now, Appendix A ensures the 
existence of a unique maximal regular solution 

u G o^((0, oo), L p (R^, R™)) n L p ((0, oo), ff 2 (R^, R™)), 

which satisfies divu = 0, thanks to the boundary condition. Splitting the spatial 
variable into a tangential part x € R™ -1 and a normal part y > as well as the 
solution as u = (v, w) into a tangential part v : [0, oo) x R" — > R™ -1 and a normal 
part w : [0, oo) x R™ — > R, we obtain by a Laplace transformation in time and a 
Fourier transformation in the tangential variable 

w 2 t;(A, £, y) - /i^fi(A, £, y) = -*£p(A, £, y), RcA > 0,£ e R"" 1 , y > 0, 

w 2 ^(A, £, !/) - K 2 *(A, C y) = -a„p(A, £, y), Re A > 0, £ G R"" 1 , y > 0, 

where we abbreviated as usual 



A £ :=e + A, u> := ^ p\ £ +mI£I 2 - 
If a = 0, then the boundary conditions read 

[v](X, 0=0, Re A > 0, £ G R™ -1 , 
^ T [w](A, + [5„i2>](A, = 0. Rc A > 0, £ € R n_1 

and we obtain 

oo 

v(X, y) = ~y G_ (A, e, y, »y) i£p(A, »?) d»j, Re A > 0, £ G R"" 1 , y > 0, 
o 

oo 

u>(A, £, y) = ~y G+(A, y, »7) 5 y p(A, rj) drj, Re A > 0, £ G R"" 1 , y > 
o 

where 

G(A, y, rj) = ^-^e~^ ly ~ vl , Re A > 0, £ G R™ -1 , y G R, »j G R 
denotes the fundamental solution to the ordinary differential equation 

^ 2 0(A, y) - /^(A, & w) = /( A > & w). Re A > o, e e R"" 1 , y e R 

and 

G±(A, y, t?) = G(A, e, y, v) ± G(A, y, -v), 

Re A > 0, £ G R n_1 , y > 0, ry > 0. 
Hence, if p is obtained as a solution to the Neumann problem (Ejv)^°' R+ , we have 

-d y p(\, y) - (-A r )- 1/2 /.(A, 0, Re A > 0, £ e R""\ y > 0, 

which implies 

oo 

N(A, = J G+(A, 0, ry) |£| e -l^ (-ArT^MA, d?7 
(To,o) o 

= Z ^ m H\Z), ReA>0, € €R-S 
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where C := y/J^^, i-e. T°[w] — h with T° as defined in Section^ If a = ±1, then 
the boundary conditions read 



T^idyv] (A, - /iV, [u)] (A, = 0, Re A > 0, £ <E 



pn— 1 



if T [u](A, + [d„it>](A, £) = 0, Re A > 0, £ e M™ _1 

and we obtain 

oo 

v(X, £, y) = -J K±(\, e, y, tj)i£0(A, £, »?)dtj, Re A > 0, £ € R™ -1 , y > 0, 
o 

oo 

w(X, e, y) = - / X*(A, £, y, v)dyP(\ 6 »?)d»?, Re A > 0, £ <G R"" 1 , y > 0, 



where 



Jf±(A, & y, »?) := f 1± — (CJ + I 9 ICI I G+(A, & y, *?) 

T 5 tr-(A, t, y, 77) 

( W + 2|C|) 2 ±|CI 2 



and 



^(A, e, * V) := 1 - ( " + 2|CI)l ^ IC ' 2 G + (A, e , y, „) 
«>v v, /; ^ (w + 2|C|) ± |Cl 2 / 

H o k-(A, C) y, 77) 

(c + 2|C|) 2 ±|C| 2 

for Re A > 0, £ € M n_1 , y > 0, i] > 0. Hence, if p is obtained as a solution to the 
Neumann problem (Ejv)^' + , we have 

-d y p(X, £, y) = (-A r y 1/2 h(\, 0, Re A > 0, £ € R"" 1 , y > 0, 

which implies 



[*&](A, = / K±(X, f, 0, 7?) |£|e-l«l" (-A r )- 1/2 /7(A, f) 
(Ti.o) o 

= ^2^ja^»0. ReA^O^eK"- 1 , 

i.e. T* 1 ^] = h with T ±:L as defined in Section H 

On the other hand, if p is obtained as a solution to the Dirichlet problem 

(Eyj)^ ' + , we have 

-d v p(X, y) = |f|e-l«l* &(A, 0, ReA > 0, C e y > 0, 
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which implies 

-2(i[d y w] + [p] 



oo 








i.e. 5 ,±1 (-2/j[a y w] + [p]) = h with S' ±1 as defined in SectionH 




n-l 
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